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Chapter  1:  Introduction 
1.0  General 

Much  interest  in  recent  years  has  centered  on  the  field  of 
adaptive  optics.  (1)  Aberrations  in  optical  systems  such  as 
telescopes  can  be  minimized  in  real  time  by  the  insertion  of  a 
phase  correcting  mirror  on  a  piezoelectric  substrate.  The 
mirror  surface  is  deformed  by  a  set  of  actuators.  Typically 
a  small  number  of  actuators  (less  than  twenty)  is  used,  each 
requiring  adjustable  voltages  on  the  order  of  several  thousand 
volts.  Such  high  voltages  are  required  to  move  the  surface 
by  O.lum  since  the  ceramic  piezoelectric  is  practically  incom¬ 
pressible. 

The  research  described  here  has  as  its  goal  the  creation 
of  adaptive  optical  components  requiring  lower  voltages,  lower 
weight,  and  lower  cost.  Also,  the  structures  under  study,  if 
successful,  would  be  capable  of  significantly  greater  phase 
control.  Thus  controllable  lenses  and  beam  formers  could  be 
created. 

A  new  area  of  polymer  technology  is  the  use  of  polymeric 
materials  in  electronic  devices  such  as  transducers,  switching 
systems,  or  electronically  tunable  lenses  (2).  Among  the  various 
polymeric  materials,  low  modulus  silicone  polymers  are  used 
for  the  majority  of  optical  devices.  The  large  deformations 
of  crosslinked  polydimethylsiloxane  (PDMS)  induced  by  uniaxial 
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compression  lead  to  applications  of  deformable  microcomponents. 

A  variety  of  optical  components  such  as  the  Ruticon  (2)  has  been 
developed  in  which  image  storage  is  obtained  by  optically 
inducing  a  charge  distribution  which  deforms  the  elastomer. 

The  metallized  elastomer  then  reflects  the  wavefront  as  shown 
in  Figure  1.0-1.  These  new  concepts  are  of  potential  importance 
in  the  field  of  imaging  and  image  display. 

1 . 1  Objectives 

One  primary  objective  of  this  study  was  the  creation  of 
a  microcomponent  consisting  of  metal  film  overlaid  with  a 
deformable  polymer  network  and  a  reflecting  metal  film  as 
shown  in  Figure  1.1-1.  Typical  dimensions  of  the  microcompo¬ 
nent  are  20  ym  in  width  and  40  ym  in  height.  The  surface 
of  the  structure  should  be  capable  of  reflecting  incident 
light  and  conducting  electricity  .  By  applying  a  voltage 
across  such  a  structure,  shown  in  Figure  1.1-1,  the  electrical 
force  induces  a  mechanical  compression  of  the  structure  due  to 
the  deformation  in  the  soft  elastomer.  The  structure  could  be 
used  as  a  focusing  mirror  or  an  electronically  controlled  beam 
scanner. 

It  is  recognized  that  the  tensile  modulus  for  this  com¬ 
pression  deformation  is  approximately  a  factor  of  three  greater 
than  that  for  a  shear  deformation.  In  order  to  take  advantage 
of  this  difference  of  moduli,  and  inverted  L-shaped  structure 
shown  in  Fig.  1.1-2  which  deforms  according  to  a  shear  mode  was 
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Fig.  1.0-1  Mechanism  of  Ruticon  Image  Formation 


When  the  photoconductor  is  illuminated,  charges 
move  toward  the  photoconduct or- elastomer  boundary. 
These  charges  modulate  the  electrical  field  and 
the  top  surface  of  the  elastomer  deforms. 
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Fig.  1.1-1  Typical  Polymeric  Microstructures 
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Side  View  of  a  Microplate  Element  Undergoing  Deflection 


Figure 
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proposed  to  maximize  the  deformation.  The  equation  of  motion 
for  the  shear  deformation  as  a  function  of  voltage  is  presented 
in  detail  in  Appendix  1A.  A  future  objective  is  to  fabricate 
such  a  structure. 

Another  area  of  research  is  concerned  with  adaptive 
optical  components  based  on  transmission  rather  than  reflection. 
A  liquid  crystal  is  capable  of  large  index  of  refraction 
changes  under  electronic  control.  We  have  studied  the  con¬ 
cept  of  a  liquid  crystal  lens  for  image  formation  as  well  as 
aberration  control. 

1.2  Photochemical  lithography  of  silicone  elastomers 

(a)  Materials  characterization 

In  this  study,  the  structure  and  properties  of  cross linked 
silicone  polymers  and  their  applications  as  microelectronic 
components  were  explored.  Equilibrium  swelling  techniques 
have  been  used  to  characterize  the  elastic  moduli  of  cross- 
linked  silicone  polymer  gels  and  the  uncrosslinked  fractions 
in  the  curing  process.  Both  large  samples  and  thin  films 
were  investigated.  The  consistency  in  uncrosslinked  fractions 
between  the  two  samples  indicated  that  both  crosslinking 
reactions  followed  the  same  mechanism  and  the  elastic  moduli 
obtained  from  large  samples  could  be  used  to  estimate  those 
of  the  microcomponents . 
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(b)  Photochemical  lithography 

A  photochemical  lithographic  technique  was  devised  for 
the  fabrication  of  a  microcomponent  on  a  glass  substrate.  This 
technique  consisted  of  the  exposure  of  the  polymer  gels  to 
ultraviolet  radiation  through  an  aluminum  mask  and  the  develop¬ 
ment  of  the  elastomer  pattern.  The  physical  shape  of  the 
developed  pattern  was  explained  in  terms  of  the  uncrosslinked 
fraction  data  obtained  from  the  equilibrium  swelling  experiments. 

(c)  Infrared  characterization 

Since  the  micro structures  are  thin  and  have  one  surface  fixed 

to  a  plate, there  is  some  difficulty  in  obtaining  equilibrium  swell¬ 
ing  measurements.  An  alternative  method  of  characterization  of 
the  crosslink  density  is  by  infrared  spectroscopy.  The  absorbance  spectra 
were  identified  according  to  the  bond  structure  of  the  silicones. 

The  bonds  reactive  in  crosslinking  were  measured  with  respect 
to  time  of  curing  to  determine  crosslinking  kinetics.  In 
addition  the  absorbance  in  the  infrared  was  found  as  a 
function  of  sample  thickness.  The  results  yielded  information 
as  to  polymer  -  polymer  interactions  and  the  effect  of  UV 
radiation  on  the  silicones. 

1.3  Metallization  of  silicone  elastomers 

For  high  quality  reproduction  of  images  these  devices 
require  large  reversible  deformations  of  the  elastomer  and  a 
highly  reflective  metallic  coating.  Ruticons  used  cross- 
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linked  silicone  elastomers  for  their  mechanical  properties 
and  various  metals  for  the  reflecting  surface.  The  amount 
of  metal  required  for  good  reflectivity  was  sufficiently  thick 
to  alter  the  mechanical  behavior  of  the  elastomer  and,  there¬ 
fore,  its  imaging  capabilities  were  limited. 

We  developed  a  cryogenic  metallization  technique  to  coat 
the  device  with  a  very  thin  reflecting  metallic  film  that  may 
not  perturb  the  elastic  properties  of  the  silicone  substrate. 

The  classical  nucleation  theory  was  used  to  interpret  this 
phenomenon . 

1.4  Theoretical  Models  of  static  deformation 

Theoretical  models  have  been  formulated  to  analyze  the  static 
deflection  of  the  polymer  microstructure.  The  equilibrium  equations 
were  derived  from  an  energy  function  composed  of  two  parts. 

One  contribution  was  the  strain  energy  as  found  from  the  theory 
of  rubber  elasticity,  the  other  that  of  the  electrostatic  energy 
due  to  a  potential  across  a  capacitor.  This  model  yielded 
expressions  for  the  static  deformation  in  terms  of  the  applied 
potential  and  the  modulus  of  the  gel. 

1 . 5  Experimental  measurements  of  the  deformation 

The  deformation  of  the  gel  was  measured  by  an  interfero¬ 
metric  technique  using  the  red  band  from  a  He  Ne  laser.  The 
structures  had  both  thin  film  and  finger  electrode  geometries 
for  exciting  the  deformation.  The  thin  film  electrodes 
yielded  interference  patterns  from  which  a  deformation  could 
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be  measured  as  a  function  of  voltage.  Both  gold  and  aluminum 
were  used  for  the  thin  film  electrodes  and  measurements  for 
the  deformation  as  a  function  of  voltage  were  obtained.  How¬ 
ever,  the  interference  pattern  for  the  finger  electrodes  was 
not  well  defined;  this  led  to  an  investigation  of  optimization 
of  parameters  in  the  metallization  technique  so  that  an  im¬ 
proved  surface  could  sustain  a  pattern. 

1.6  The  dynamic  characterization  of  polymer  films 

In  order  to  model  the  dynamics  of  the  structure  and  esti¬ 
mate  response  times,  the  viscous  and  elastic  tensor  components 
have  to  be  determined.  Since  the  gel  is  isotropic,  there  are 
only  four  constants  to  be  determined;  two  elastic  and  two  viscous. 
A  non- destructive  experimental  technique  to  determine  these 
constants  is  the  propagation  of  plate  modes  in  a  delay  line 
apparatus.  The  dispersion  relations  for  lossy  plate  modes  were 
derived  for  an  anisotropic  film  of  which  the  gel  is  a  special 
case.  An  example  of  the  evaluation  of  elastic  and  loss  constants 
was  made  for  PVF2. 

1.7  Liquid  crystal  lens  -  Analysis 

We  have  analysed  the  concept  of  an  electronically  controlled 
lens  based  on  the  concept  of  index  of  refraction  modulation  by 
an  electric  field.  The  materials  required  and  the  possible 
geometries  of  the  structure  have  been  studied.  The  optical 
performance  has  been  evaluated  assuming  simple  models  for  the 
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crosstalk  between  electrodes.  For  lens  applications,  this 
crosstalk  (coupling)  greatly  improves  the  optical  performance 
by  extending  the  spatial  frequency  response. 

1.8  Liquid  crystal  lens  -  Experimental  effort 

We  have  begun  the  experimental  phase  by  procuring  a  com¬ 
mercial  liquid  crystal  cell  and  modifying  it  for  our  purposes. 

No  conclusive  results  have  been  obtained  so  far. 

1.9  Pattern  generator 

We  have  nearly  finished  the  construction  of  our  pattern 
generator.  A  5-watt  Spectra  Physics  argon  ion  laser  is  the 
optical  source  for  the  system  which  is  composed  of  three  (x,  y, 

0)  stepper  motor-driven  platforms  under  microcomputer  control. 
This  system  is  capable  of  1  ym  precision  with  the  potential  of 
a  1  ym  spot  size.  The  laser  can  be  operated  in  the  visible 
(green,  blue)  region,  or  in  the  near  ultraviolet.  The  system 
is  operational  and  is  being  calibrated  and  adjusted  to  maximize 
accuracy  and  minimize  vibration.  The  computer  has  been  pro¬ 
grammed  for  simple  patterns  such  as  straight  lines. 

This  system  will  permit  selective  crosslinking  and  pyrolysis 


of  the  elastomer  materials. 


1.10  Statement  of  Goals 


The  specific  goals  of  the  research  were: 

(1)  to  investigate  the  preparation  and  characterization  of 
polymer  gels  for  electronic  applications  (Sections 
2.1,  2.2,  and  2.4); 

(2)  to  develop  techniques  for  the  lithography  and  metalli¬ 
zation  of  these  materials  (Section  2.3); 

(3)  to  determine  the  effects  of  an  electric  field  on  the 
viscoelastic  behavior  of  the  components  (Sections  2.5 
and  2.6) ; 

(4)  to  develop  a  technique  to  find  the  viscoelastic  constants 
to  be  used  for  modelling  a  general  structure  (Section  3.0); 

(5)  to  analyze  and  fabricate  an  adaptive  electro -optic  lens 
based  on  liquid  crystal  technology  (Section  4) . 
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Chapter  2  Elastomer  Structures 
2.1  Materials 

The  materials  selected  for  this  study  were  low  modulus  silicone 
gels.  These  materials  were  chosen  because  of  their  low  modulus  and 
their  dimensional  and  thermal  stability.  Organic  elastomers  such 
as  poly (isoprene)  do  not  have  the  required  thermal  stability  and 
very  low  modulus  necessary  for  these  applications. 

The  silicone  is  to  be  used  as  a  gel,  that  is,  the  individual 
molecules  are  chemically  bonded  to  each  other  until  the  gel  point 
is  attained, i.e. ,  a  single  molecule  is  formed.  The  mechanisms 
and  techniques  of  cross linking  and  characterization  of  the  net¬ 
works  are  presented  in  this  section. 

A  silicone  polymer  is  composed  of  the  repeating  unit: 


and  R2  are  commonly  chosen  to  be  one  of  the  following  groups, 
CH2  -  CH,  CH3,  CF3CH2CH2,  CNCH2CH2.  For  Rx  -  R2  -  CH3,  poly- 
dimethylsiloxane  results.  The  crosslinking  begins  by  initiating 
a  radical  (unbonded  electron)  on  the  methyl  group.  This  occurs 
when,  upon  heating,  a  peroxide  initiator  (mixed  into  the  polymer) 
breaks  into  a  diradical  by  the  reaction: 
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14 0°C  . 

R  -  R  +  Q  -►  R+R 


The  radical  peroxide  transfers  the  electron  to  a  methyl  site: 


The  unpaired  electron  of  the  side  group  carbon  atom  is  now 
reactive  and  capable  of  bonding  to  a  similar  species  as  follows 


The  resulting  crosslink  is  bracketed  above. 

The  crosslinking  occurs  randomly  in  the  network;  however, 
fluctuations  may  occur  due  to  aggregation  of  the  peroxide  so 
that  a  nonuniform  crosslinking  density  occurs.  This  problem 
may  be  circumvented  by  periodically  substituting  a  diene  CH2  “ 
CH  group  for  a  methyl  group  in  the  polymer  chain.  Sine  the 
diene  group  is  more  susceptibel  to  radical  formation,  uniformly 
distributed  diene  groups  result  in  a  network  with  a  uniform 
crosslink  density.  For  the  vinyl  crosslinking  an  initiator 
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which  is  unreactive  v?ith  saturated  bonds  and  reactive  only  with 
the  double-bonded  vinyl  group,  such  as  dicumyl  peroxide,  is  used. 

2 . 2  Characterization  of  the  Network  or  Modulus  Determination 

The  density  of  crosslinks  or,  equivalently,  the  number  of 
monomer  units  per  network  chain  may  be  found  by  equilibrium 
swelling  and  by  stress-strain  measurements.  For  small  strains, 
the  tensile  modulus,  E,  of  the  network  is 

E  =  3RTve  (2.2-1) 

where  R  is  the  gas  constant,  T  is  the  absolute  temperature,  and 
vg  is  the  number  of  chains  per  unit  volume.  The  density  of  cross¬ 
links  may  be  determined  using  this  linear  relation.  Equation 
2,2-1  for  the  modulus  was  derived  from  a  kinetic  theory  which 
assumes  no  interaction  between  chains.  It  is  known  that  the  re¬ 
tractive  stress  possesses  an  additional  term  which  has  been 
physically  assigned  to  entanglements  and  transients  of  polymer 
chains . 

This  entanglement  term  has  been  shown  to  go  to  zero  when 
the  network  is  swollen.  We,  therefore,  outline  the  characteriza¬ 
tion  of  silicone  by  swelling.  A  detailed  description  of  the 
experimental  procedure  in  the  swelling  experiment  is  given  in 
Appendix  2A. 

The  numerics  of  the  cross  linking  density  are  given  in  this 
section.  We  wish  to  calculate  the  number  of  chains  in  a  network. 
Consider  a  swollen  silicone  gel  in  equilibrium  with  a  solvent 
(hexane).  The  volume  fraction  of  the  polymer,  V£.  is 
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v2  -  -L  /  C(w/ph)  +  (l/pr)  (2.2-2) 

where  pr  is  the  density  of  the  polymer,  p^  is  the  density  of 
the  solvent  (hexane) ,  and  w  is  the  grams  of  hexane  absorbed 
per  gram  of  polymer.  The  number  of  modes  of  effective  network 
chains  per  cubic  centimeter,  vg,  for  tetrafunctional  sites  is 
given  by  the  relation  ^ 


C1n(1“V  +  V2  +  X1V2] 


V, 


(vF^  V 


2/2 


) 


(2.2-3) 


where  is  the  solvent /polymer  interaction  parameter  (0.48  for 
hexane) ^  and  is  the  molar  volume  of  haxane  (130.6  cc/mole) . 
The  value  of  vg  can  be  compared  to  the  modulus  since,  from  Eq. 
2.2-1. 

E 

ve  3RT 

This  determination  of  the  modulus  depends  only  upon  the  number 
of  true  network  chains  and  is  not  affected  by  entanglements. 
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2.3  Fabrication  of  Structures 


The  preparation  and  the  lithography  of  elastomers  based 
on  polydimethylsiloxane  were  investigated.  The  silicone 
materials,  crosslinked  either  by  thermal  or  ultraviolet 
radiation,  were  characterized  by  equilibrium  swelling 
measurements . 

A  process  was  devised  for  the  fabrication  of  an  elastomer 
microstructure  on  a  glass  substrate  through  the  use  of  semi¬ 
conductor  photolithographic  techniques.  This  technique  con¬ 
sisted  of  the  exposure  of  the  polymer  to  ultraviolet  radiation 
through  an  opaque  mask,  the  development  of  an  elastomer  pattern, 
and  the  metallization  of  the  polymer  surface. 

A  microcomponent  was  constructed  consisting  of  a  metal 
film  overlaid  with  a  deformable  polymer  network  and  a  reflecting 
metal  film  as  shown  in  Figure  1.1-1,  Typical  dimensions  of  the 
microcomponent  are  20  urn  in  width  and  40  ym  in  height.  The 
surface  of  the  structure  should  be  capable  of  reflecting  in¬ 
cident  light  and  conducting  electricity. 

To  this  end,  a  novel  cryogenic  metallization  technique 
was  introduced  to  coat  the  surface  of  the  device  with  a  thin 
reflective  metal  layer.  The  conductance  of  the  metal  film 
was  used  to  characterize  the  metallized  surface. 

2.3.1  Materials  for  Microstructures 

The  three  materials  selected  for  this  study  were  two  experi¬ 
mental  materials  from  Dow  Corning  Corporation,  Midland,  Michigan, 
designated  by  the  supplier  as  System  A  and  System  B,  and  the 
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third  a  Dow  Corning  Dielectric  Gel.  System  A  and  System  B  are 
ultraviolet  curable  while  the  Dielectric  Gel  is  thermally 
curable.  All  three  systems  are  essentially  polydimethyl- 
siloxane  (PDMS)  with  different  functional  groups  for  cross- 
linking.  System  A  consists  of  a  mercapto-functional  siloxane 
(X2  -  7077)  and  a  vinyl-functional  fluid  (X2  -  7078)  with  the 
crosslinking  reaction  between  these  two  sites.  System  B 
is  a  vinyl-functional  siloxane  (X2  -  7079)  and  a  silane- 
functional  polymer  (X2  -  7080) .  The  Dielectric  Gel  (Q^  -  6537) 
is  a  two-part  system  consisting  of  a  vinyl  fluid  and  a 
silane  fluid.  The  mixing  specifications  for  the  components  are 
presented  in  Appendix  2B. 

2.3.2  Characterization  by  Swelling  of  Bulk  and  Microstructures 

(a)  Large  Samples  (bulk  swelling) 

The  ultraviolet  curable  samples  (Systems  A  and  B)  were 
prepared  according  to  the  supplier's  instructions  (3)  and 
exposed  to  a  200  watt  Hg  arc  light.  The  exposure  energy  of 
the  ultraviolet  light  was  measured  with  a  Laser  Precision  Corp. 
KT-4000  radiation  detector.  For  the  Dielectric  Gel,  the  samples 
were  prepared  in  aluminum  cups  and  cured  in  an  oven  at  413  K. 
After  preparation,  the  samples  were  swollen  to  equilibrium 
for  two  days, removed  from  hexane,  and  weighed.  Vaporization 
of  hexane  during  the  weighing  operation  was  minimized  by 
saturating  the  balance  atmosphere  with  hexane  vapor  before 
weighing.  The  samples  were  then  dried  in  a  vacuum  oven  for 
24  hours  at  ambient  temperature  and  reweighed  in  order  to 
determine  the  soluble  fraction  of  fluid. 
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(b)  Thin  Films  (microswelling) 

The  ultraviolet  curable  Systems  A  and  B  were  spin-coated 
on  a  glass  substrate  and  exposed  to  ultraviolet  radiation  for 
specified  periods  of  time.  The  samples  were  then  swollen  in 
an  isooctane  bath  for  one  minute,  removed,  dried,  and  weighed. 

The  same  procedure  was  repeated  three  times  for  all  samples . 

The  soluble  fraction  or  the  percentage  extractable  was  obtained 
by  measuring  the  weight  difference  before  and  after  swelling. 

The  swelling  data  were  converted  into  crosslink  densities 
and  elastic  moduli  using  the  theory  of  rubber  elasticity  (Sec. 
2.2).  The  soluble  fractions,  the  crosslink  densities,  and  the 
elastic  moduli  of  the  three  systems  are  listed  in  Tables  2.3-1, 

II,  and  III.  The  results  are  plotted  as  a  function  of  exposure 
energy  (or  curing  time)  in  Figures  2.3-1,  2.3-2,  2.3-3,  and  2.3-4. 

(c)  Soluble  Fraction  in  Large  Samples  (bulk  swelling) 

For  System  A,  the  soluble  fractions  or  the  percentage  ex- 
tractables  of  all  samples  were  157.  ±  27.  for  curing  times  greater 
than  2  minutes.  This  indicates  that  a  completely  crosslinked 
network  is  already  formed  after  exposure  to  ultraviolet  radiation 
for  only  2  minutes.  For  System  B  and  the  Dielectric  Gel,  the 
crosslink  reaction  was  not  complete  since  the  soluble  fraction 
decreased  with  increasing  curing  time.  Similar  results  were 
obtained  by  Schwartz  et  al.  (4)  on  other  silicone  polymers.  The 
scatter  which  occurred  in  the  data  for  System  B  may  be  due  to 
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Table  2.3-1 


Curing  Time (Min) 

Bulk  Swelling 

7o  Extractable 

Data  for  System  A 

v  (10^k  mole/m'*) 
e  g  .  _ 

E(MPa) 

2 

15.1 

1.49 

1.1 

3 

12.2 

1.42 

1.1 

4 

10.6 

1.78 

1.3 

6 

13.8 

1.71 

1.3 

7 

11.5 

1.83 

1.4 

8 

12.8 

1.59 

1.9 

8 
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Table  2.3-II 


Curing  Time (Min) 

Bulk  Swelling  Data  for  System 

3 

7.  Extractable  v  (k  mole/m  ) 

® _ £ _ _ 

8 

39.4 

12.5 

10 

43.3 

18.0 

12 

26.3 

18.2 

14 

28.4 

19.0 

16 

23.1 

32.5 

18 

17.1 

47.6 

20 

20.4 

11.4 

25 

22.5 

62.4 

30 

21.4 

38.6 

35 

18.9 

52.9 

40 

17.8 

10.5 

60 

14.0 

87.6 

E(Pa) 

8.3  x  104 

1.3  x  105 

1.4  x  105 

1.4  x  105 

2.4  x  105 

3.5  x  105 

8.5  x  105 
4.8  x  105 

2.7  x  105 
4.0  x  105 

7.8  x  105 

6.6  x  105 
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Table  2.3-III 

Bulk  Swelling  Data  for  Dielectric  Gel 


Curing  Time (min) 

10 

15 

20 

25 

30 


%Extractable  u 

-  — 6- 

48.5 

43,3 

39,2 

36,7 

32 

22,9 


iole/m^) 

E(104Pa) 

4.0 

2.9 

6,7 

5.0 

8,1 

6.1 

6.5 

4.9 

11.7 

8.8 

24.8 

19.0 

35 


EXPOSURE  ENERGY (Joules/cm2) 


Extractable  versus  curing  time  (or  exposure 


DIELECTRIC  GEL(I40°C) 
SYSTEM  B  (U.V.) 
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Figure  2.3-2#  Extractable  versus  curing  time  (or  exposure 
energy)  for  System  B  and  Dielectric  Gel 


EXPOSURE  ENERGY  ( Joules/cm 


%  Extractable  versus  curing  time  (or  exposure 
energy  )  for  System  B  • 
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U.V  EXPOSURE  ENERGY  (Joules/cm2 ) 


Figure  2,3-4  Variation  of. elastic  modulus  with  curing 
time  (or  exposure  energy) 
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(i)  the  variations  in  the  sample  geometry,  which  could  be  as 
high  as  5%; 

(ii)  uncertainties  in  the  swelling  measurements; 

(iii)  the  uncertainties  in  the  elastic  analysis,  which  have  been 
attributed  to  the  neglect  of  two  terms.  One  is  a  factor 
relating  the  modulus  to  an  ilolated  unperturbed  polymer 
chain.  The  second  effect  is  that  of  chains  which  are 
ineffectively  crosslinked  (i.e.,  crosslinked  at  only 

one  end)  and  reduce  the  modulus . 

(d)  Crosslink  Densities  and  Elastic  Moduli  of  Large  Samples 

As  shown  in  Figure  2.3-4,  the  elastic  modulus  of  System  A 
was  almost  constant  with  curing  time.  This  consistent  with  the 
soluble  fraction  data  since  System  A  formed  a  completely  cross¬ 
linked  network  after  an  exposure  of  only  2  minutes  whereas ,  for 
System  B  and  the  Dielectric  Gel,  the  crosslinking  reaction  was 
not  complete  since  the  elastic  modulus  increased  with  increasing 
curing  times  which  afrees  with  the  results  of  Schwartz  et  al .  (4) . 

(e)  Thin  Films  (microswelling) 

The  objective  of  microswelling  was  to  determine  the  uncross- 
linked  soluble  fraction  in  the  microstructure.  The  results 
were  then  used  to  estimate  the  modulus  of  the  microelectronic 
components  from  the  microstructure  thickness  after  developing. 

For  the  ultraviolet-curable  Systems  A  and  B,  the  soluble  fractions 


are  listed  in  Tables  2.3-IV  and  2.3-V,  and  the  results  are  plotted 
in  Figures  2.3-1  and  2.3-3.  For  both  systems,  the  samples  reached 
a  constant  weight  after  immersion  for  three  minutes  in  isooctane 
indicating  that  this  time  was  sufficient  for  completely  developing 
the  microstructures.  As  shown  in  Figure  2.3-3,  the  percent  extract 
ables  for  System  B  obtained  from  microswelling  were  consistent 
with  those  from  bulk  swelling.  Hence,  the  elastic  moduli  obtained 
from  large  samples  could  be  used  to  estimate  the  moduli  of  the 
microstructures . 

2.3.3  Photolithography  and  Sample  Preparation 

As  stated  in  Section  2.3  the  specific  objective  of  the 
study  was  the  creation  of  a  micro component  consisting  of  a 
metal  film  overlaid  with  a  polymer  network  and  a  reflecting  metal 
film  as  shown  in  Figure  1.1-1.  When  an  external  voltage  is 
applied  across  the  structure  in  the  deformation  testing,  the  ex¬ 
ternal  voltage  will  be  shorted  if  metal  is  deposited  on  the  ver¬ 
tical  edges.  Only  a  vertical  or  an  undercut  structure  can  be 
metallized  without  depositing  metal  vapor  on  the  vertical  edges. 

For  this  purpose,  a  photochemical  lithography  process  was  devised 
for  the  fabrication  of  microcomponents  on  a  glass  substrate  (5,6). 
The  process  included  spin-coating,  exposure,  and  developing.  The 
physical  shape  of  the  structure  was  explained  in  terms  of  the 
uncrosslinked  fraction  data  obtained  from  the  equilibrium  swelling 
experiments . 

The  photolithographic  process  first  involves  pretreatment 
of  the  glass  substrate.  Next,  a  photosensitive  polymer  film  (photo 
resist)  is  applied  to  the  substrate  by  spin-coating  techniques. 
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Table  2.3-IV 

Microswelling  Data  for  System  A 


Curing  Time 

Weight  of  Film(g) 

%Extractable 

Min. 

0  min, 

.  1  min. 

2  min. 

3  min. 

2 

.0339 

.0243 

.0243 

.0218 

35.6 

4 

.0335 

.0280 

.0274 

.0274 

18.2 

5 

.0345 

.00317 

.0305 

.0305 

11.6 

6 

.0361 

,0331 

.0331 

.0306 

15 

7 

.0350 

.0313 

.0316 

.0313 

10.6 

8 

.0370 

,0348 

.0336 

,0336 

9.1 

3 

* 

•¥ 

5 


i 

I 

1 


30 


Table  2.3-V 

Microswelling  Data  for  System  B 

Curing  time (Min)  Weight  of  Film(g)  %  Extractable 


min 

0  min 

1  min 

2  min 

3  min 

8 

.0281 

.0176 

.0165 

.0165 

41.3 

10 

.0269 

.0179 

.0179 

.0179 

33.5 

12 

.0266 

.0196 

.0178 

.0178 

33 

16 

.0271 

.0211 

.0201 

.0201 

26 

20 

.0289 

.0231 

.0223 

.0223 

22.8 

25 

.0281 

.0235 

.0230 

.0230 

18.2 

30 

.0282 

.0239 

.0239 

.0239 

15.3 

35 

.0252 

.0218 

.0218 

.0218 

13.7 

40 

.0269 

.0236 

.0236 

.0227 

15.5 

60 

.0270 

.0245 

.0240 

,0236 

12.5 

60 


0270 


0245 
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The  film  is  then  selectively  irradiated  with  an  ultraviolet  lamp 
through  a  mask  which  contains  clear  and  opaque  features  outlining 
the  desired  pattern.  The  areas  in  the  photoresist  that  are 
exposed  to  light  are  made  either  soluble  or  insoluble  to  a  solvent 
(developer)  (6).  For  a  positive  photoresist , the  irradiated 
regions  dissolve  leaving  a  positive  image  of  the  mask  on  the 
substrate  while,  for  a  negative  resist,  the  unexposed  regions 
dissolve  leaving  a  negative  image  of  the  mask  on  the  substrate. 

The  photochemical  lithography  process  used  in  this  study 
is  diagrammed  in  Table  2,1, 
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2.3*4  Characterization 

A  conventianal  polymer  resist  pattern  was  prepared  and 
metallized  as  shown  in  Figure  2.3-5.  The  resist  pattern  was 
then  removed  by  acetone  and  the  resulting  aluminum  pattern  was 
used  as  a  mask  for  the  polymer  gel  pattern.  For  the  deformable 
microstructure  fabrication,  a  very  thin  transparent  aluminum 
film  (6  nm)  was  evaporated  over  the  opaque  pattern  as  shown 
in  Figure  2.3-6.  The  polymer  gels  were  then  spun  over  the 
aluminum  pattern  as  shown  in  Figure  2.3-7.  The  exposed 
area  crosslinked  while  the  unexposed  area  remained  uncross- 
linked.  Since  the  crosslinked  gel  is  insoluble,  the  soluble 
fraction  was  washed  away  by  a  solvent. 

The  relation  between  the  microcomponent  thickness  and 
spinning  speed  was  investigated  by  spin-coating  the  polymer 
gels  en  glass  substrates  at  various  speeds.  The  spin-coating 
results  are  given  in  Appendix  2C.  The  polymer  films  were  thin 
cured  in  a  200  watt  mask  aligner.  The  thickness  was  obtained 
by  measuring  the  cross  section  with  an  optical  microscope. 

The  microelectronic  patterns  were  investigated  by  electron 
and  optical  microscopy.  To  facilitate  the  electron  microscope 
investigation,  a  thin  layer  of  aluminum  was  evaporated  upon  the 
sample  under  high  vacuum  conditions  to  prevent  charge  build-up 
on  the  polymer  surfaces. 

Scanning  electron  micrographs  of  the  electronic  components 
prepared  with  System  B  are  shown  for  the  20  pm  wide  pattern  in 


METALLIZED  POSITIVE  PHOTORESIST  PATTERN 

UZ7Z&  ALUMINUM 
I  I  PHOTORESIST 


Fig.  2.3-5  Photoresist  Pattern  for  Lift  Off  Technique 


ALUMINUM  MASK 


YZTZtt  OPAQUE  ALUMINUM  PATTERN 
TRANSPARENT  ALUMNUM 


Fig.  2.3-6  Aluminum  Contact  and  Mask  Layer  For 
Microstructure  Fabrication 
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PHOTOCHEMICAL  LITHOGRAPHY 


i 

i 

u.v. 

0  .*  Crosslinked  Silicone 

\ 

□  *.  Uncrosslinked  Silicone 


Fig.  2.3-7  Microstructure  Photocythography  Process 
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Figure  2.3-6  and  for  the  60  urn  wide  spacing  in  Figure  2.3-9. 

The  polymer  thickness  in  both  samples  is  13  pm.  Due  to  the 
lithography  process,  the  narrow  pattern  is  more  rectangular 
than  the  wide  pattern.  in  developing  the  pattern,  the  solvent 
extracts  the  unexposed  area  and  the  extractable  fraction  of 
the  exposed  area.  The  exposed  area  contracts  after  drying 
but  this  contraction  can  occur  on  only  three  of  the  four 
sides  of  the  pattern  since  the  adhesion  between  the  polymer 
and  the  substrate  is  strong  enough  to  hold  the  sample  onto 
the  substrate.  Hence,  by  determining  the  cross  sectional 
areas  of  the  pattern  and  a  non-extracted  pattern,  the  extract- 
able  fraction  of  the  microstructure  can  be  estimated  and  compated 
to  the  equilibrium  swelling  results.  This  is  shown  in  Figure 
2.3-10  in  which  an  undeveloped  and  a  developed  pattern  are  shown. 
For  this  sample,  an  extractable  fraction  of  237.  was  obtained 
which  is  consistent  with  the  equilibrium  swelling  measurements 
on  bulk  samples  at  similar  curing  times. 

An  optical  micrograph  of  the  electronic  component  prepared 
with  System  B  is  shown  in  Figure  2.3-H*  The  polymer  was 
cured  for  60  minutes  and  has  a  thickness  of  35  pm.  This 
pattern  is  much  more  rectangular  than  the  13  pm  thick  pattern. 

The  structure  can  be  used  to  construct  a  microelectronic  compo¬ 
nent  as  shown  in  Figure  1.1-1. 


2.3.  5  Metallization  of  Polymer  Gels  in  the  Rubbery  State  and 
the  Glassy  State 

It  was  pointed  out  in  the  previous  chapter  that  the  micro¬ 
structure  will  be  used  as  a  fecusing  mirror  or  an  electronically 
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Figure  2.3-9 


Elastomer  pattern  with  60  ym  spacing. 


I 


Fig.  2 
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Figure  2.3-H  Elastomer  pattern  for  deformation  testing. 
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controlled  beam  scanner  (7) .  The  surface  of  the  structure  should 
be  capable  of  reflecting  incident  light  and  conducting  electricity. 
For  these  purposes,  a  reflective  and  conductive  thin  metal  film 
was  coated  on  the  top  surface  of  the  components.  The  metal  films 
deposited  at  ambient  conditions  with  high  reflectivity  were  thick 
enough  to  alter  the  mechanical  behavior  of  the  elastomer  and, 
therefore,  its  imaging  capabilities  were  very  limited.  A  cryogenic 
metallization  technique  was  developed  to  coat  the  component 
surfaces  with  a  thin  reflective  metal  layer  (8) .  The  metallization 
process  was  interpreted  in  terms  of  the  physical  properties  of 
polymer  gels  and  the  classical  thin  film  nucleation  theory. 


2.3*6  Experimental  Procedures  for  Metallization 

The  materials  were  mixed  according  to  specification  and 
spin^cast  on  glass  microscope  slides.  Samples  of  four  different 
crosslink  densities  were  prepared  for  each  system.  For  System  A, 
the  samples  were  exposed  to  a  200  watt  mask  aligner  for  0,  1,  2, 
and  4  minutes  while,  for  System  B,  the  samples  were  exposed  for 
0,  5,  10,  and  15  minutes.  Shorter  exposure  times  for  System  A 
were  chosen  since  the  kinetics  of  crosslinking  was  faster  than 
for  System  B.  For  the  Dielectric  Gel,  one  sample  was  cured  at 
room  temperature  for  25  hours  and  the  others  were  cured  in  an 
oven  at  413  K  for  10,  20,  and  30  minutes. 

For  the  purpose  of  cooling  the  polymer  gels  below  their 
glass  transition  temperatures,  a  low  temperature  metallization 
vacuum  chamber  was  fabricated  as  shown  in  Figure  2.3-12.  The 
apparatus  was  constructed  of  a  glass  bell  jar  with  a  two-piece 
metal  cover.  The  top  consisted  of  an  aluminum  plate  fitted  with 
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Fig. 


THERMOCOUPLE  WIRE 


BASE  PLATE 


2.3-12  Vacuum  Chamber  for  Low  Temperature  Metal  Evaporation 
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a  hollow  cylindrical  stainless  steel  tube  with  a  brass  cover. 


When  liquid  nitrogen  is  placed  in  it,  the  tank  acts  as  a  low 
temperature  reservior  for  the  polymer  samples  clamped  in  thermal 
contact  to  it. 

For  every  evaporation,  four  samples,  each  of  a  different 
crosslink  density,  were  simultaneously  metallized  to  maintain 
a  uniform  thickness  from  sample  to  sample.  It  was  found  for 
both  indium  and  altiminum  that  a  0.635  cm  length  of  0.127  cm 
diameter  wire  evaporated  at  2  x  10”^  torr  was  sufficient  for  the 
formation  of  a  metal  film.  The  films  were  about  50  pm  in  thick¬ 
ness.  For  the  low  temperature  evaporation,  the  silicone  was 
allowed  to  equilibrate  for  30  minutes  at  low  temperature. 

After  ten  minutes  of  equilibration  time,  the  temperature  of 
the  silicone  was  measured  to  be  100  K  using  a  thermocouple  in¬ 
serted  into  the  vacuum  chamber,  thus  assuring  that  the  polymer 
was  well  below  its  glass  transition  temperature  of  approximately 
150  K.  (9)  • 

The  resulting  metal  coatings  were  investigated  by  optical 
microscopy  and  measurements  of  reflectance  and  conductance.  The 
conductance  was  measured  by  a  two-point  microprobe  with  a  separa¬ 
tion  distance  of  0.6  cm.  For  the  reflectance  measurement,  a 
He  Ne  laser  beam  was  used  to  scan  a  portion  of  the  metal  film 
that  had  the  highest  reflectance.  The  selected  area  was  scanned 
in  a  horizontal  motion  for  three  vertical  heights  and  the 
largest  value  of  the  ref lectance >as  measured  by  a  CdS  cell,  was 
recorded.  The  sample  averages  were  normalized  to  the  reflectance 
measured  from  an  aluminum  metallized  glass  slide. 


■*> VV::  ~ 


2.3.7  The  Characterization  of  the  Aluminum  Metallization 
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In  this  and  the  following  section,  the  results  Of  the  metal 
surface  morphology  for  aluminum  for  indium  are  presented.  We  com¬ 
pare  the  surfaces  for  both  low  temperature  and  ambient  temperature 
thermal  processing. 

The  reflectivities  and  the  conductances  are  listed  in  Tables 
2.3-VI  and  2.3-VII.  All  of  the  aluminum  surfaces  formed  by  evapor¬ 
ation  onto  the  substrate  at  ambient  temperature  appeared  opaque  with 
an  oxidized  surface.  Optical  reflective  microscopy  showed  the 
surface  texture  to  consist  of  a  fine  structure  of  folds  with  no 
long-range  correlation.  This  metallic  texture  is  shown  in  micro¬ 
graph  (A)  of  Figure  2.3-13. 

As  a  result  of  the  low  temperature  evaporation,  the  gels  of 
all  three  systems  had  a  mirror- like  surface.  Optical  microscopy 
for  these  samples,  for  a  lower  magnification  than  those  of 
ambient,  shown  in  Figure  2.3-13-B,  revealed  an  aluminum  film  that 
contained  wrinkles  with  a  longer  range  of  order  as  compared 
to  samples  produced  at  ambient  temperature.  Both  the  aluminum 
and  the  indium  surfaces  exhibited  some  features  in  common  for  the 
low  temperature  processing.  First,  large  wrinkles  which  may  be 
caused  by  flow  could  be  seen  for  uncured  samples  (unexposed  to 
ultraviolet  radiation).  Second,  microscopy  of  the  mirror-like 
surfaces  of  the  low  temperature  gels  showed  that  a  less  coarse 
wrinkled  texture  occurred  for  the  cured  materials.  The  size  of 
the  wrinkles  decreased  with  increasing  crosslink  density  of  the 
gel  as  shown  in  Figure  2.3-14.  Unlike  the  ambient  samples,  the 
surfaces  for  these  samples  had  cracks  which  are  probably  due 
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A.  Aluminum  evaporated  onto  an  ambient  temperature  substrate. 
The  gel  is  system  A  cured  for  2  min.;  the  magnification  is 
50  x  10. 


B.  Aluminum  evaporated  onto  a  liquid  nitrogen  cooled  substrate. 

The  gel  is  the  same  system  and  cure  as  above;  the  magnification 
is  10  x  10. 


Figure  2.3-13  A  comparison  of  A1  surfaces  for  low  and 
ambient  temperature  of  the  substrate 
during  evaporation. 
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Table  2.3-VI 

Surface  Characterization  of  Silicone  Gels  with  Aluminum  Metallization 


Polymer  System 

Cure  Time 
(min.) 

Normalized  Reflectance  (a) 
Metallization  Temp. 

300  K  below  T 

_ g _ _ 

_2 

Conductance (10  s  )  (b) 
Metallization  Temp. 

300  K  below  T 

R 

A 

1 

0 

99 

3 

12 

2 

61 

98 

3 

18 

4 

107 

104 

6 

20 

B 

5 

93 

99 

3 

10 

10 

32 

101 

5 

10 

15 

50 

103 

10 

15 

Dielectric  Gel 

10 

— 

— 

0 

32 

20 

17 

21 

0 

30 

30 

87 

75 

0 

37 

*  These  reflectances  were  normalized  to  that  of  a  metallized  glass  slide. 

a)  For  those  normalized  reflectance  values  within  25%  of  that  of  glass,  specular 
reflection  occurred. 

—8 

b)  No  conductance  occurred  (less  than  10  S)  for  ambient  aluminum  metallization 
onto  the  dielectric  gel.  Although  there  was  spread  in  the  values,  the  order 
of  magnitude  of  the  values  presented  are  reliable. 


48 


Table  2.3-VII 


Surface 

Characterization  of 

Silicone 

Gels  with  Indium  Metallization 

Polymer  System 

Cure  Time 
(min.) 

* 

Reflectances  (a) 
Metallization  Temp. 

300  K  below  T 

8 

—2 

Conductance (10  S  )  (b) 
Metallization  Temp. 

300  K  below  T 

8 

A 

1 

— 

— 

0 

26 

2 

112 

66 

0 

38 

4 

99 

95 

0 

43 

B 

5 

33 

87 

0 

12 

10 

70 

100 

0 

15 

15 

97 

104 

0 

20 

Dielectric  Gel 

10 

94 

90 

0 

32 

20 

94 

76 

0 

30 

30 

91 

24 

0 

37 

*  These  reflectances  were  normalized  to  that  of  a  metallized  glass  slide. 

a)  For  those  normalized  reflectance  values  within  25Z  of  that  of  glass,  specular 
reflection  occurred. 

-8 

b)  No  conductance  occurred  (less  than  10  S)  for  ambient  Indium  metallization  on 
any  of  the  polymer  substrates.  Akthough  there  was  spread  in  the  values  the 
order  of  magnitude  values  presented  are  reliable. 


t 

1 

t 
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to  the  difference  in  expansivity  between  the  polymer  and  the 
aluminum  as  the  sample  is  taken  from  liquid  nitrogen  to  room 
temperature . 

The  reflectances  listed  in  Tables  2.3-VI  and  2.3-VII  show 
some  variation.  The  higher  values  were  obtained  when  the  reflected 
laser  beam  could  be  focused  onto  the  detector.  The  lower  reflec¬ 
tance  values  occurred  when  the  reflected  beam  was  diffused,  thus 
indicating  the  smoothness  of  the  texture  of  the  metallic  surface. 

A  comparison  of  the  reflectance  for  System  A  illustrates  that  a 
higher  quality  reflecting  surface  corresponds  to  the  long-range 
correllated  texture  for  the  low  temperature  process . 

Conductances  for  all  three  systems  of  the  gels  (Tables  2.3-VI 
and  2.3-VII)  were  greater  for  the  liquid  nitrogen  cooled  gels  than 
those  of  the  ambient.  For  both  of  the  ultraviolet-cured  systems, 
the  ambient  conductances  are  a  small  fraction  of  those  prepared 
at  the  lower  temperature.  The  Dielectric  Gel  processed  at  low 
temperature  conducted  for  both  metals ,  as  contrasted  with  that 
prepared  by  the  ambient  evaporation  where  no  conduction  was  measured. 
This  is  a  striking  result  when  compared  to  the  ultraviolet  cured 
materials . 

2.3.8  The  Characterization  of  the  Indium  Metallization 

The  indium  deposition  for  ambient  conditions  yielded  a  brown- 
colored  film  which  was  transparent.  For  System  A,  with  curing 
times  of  1  and  2  minutes,  a  variety  of  colors,  which  ranged  from 
white  to  pink  to  purple  were  observed  on  the  surface.  For  low 
temperature  evaporation,  the  indium  had  a  mirror-like  surface 
and,  unlike  aluminum,  did  not  appear  to  have  any  cracks. 
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The  reflectances  of  the  indium  surface  (Table  2.3-VII)  were 
high  for  both  processing  temperatures  for  System  A,  while  for 
System  B,  the  low  temperature  process  yielded  a  greater  average 
reflectance. 

The  results  for  the  conductance  for  indium  were  similar  to 
those  of  aluminum  in  that  none  of  the  gels  conducted  under  the 
condition  of  ambient  evaporation,  whereas,  for  the  low  tempera¬ 
ture  evaporations,  the  metal  films  were  conducting. 

Little  variation  was  noted  in  the  conductivity  for  dif¬ 
ferent  crosslink  densities.  It  appears  that  there  may  be  a 
greater  surface  migration  for  indium  than  for  aluminum  and 
that  its  interaction  with  silicone  allows  a  similar  film  to  be 
formed  on  all  of  the  silicone  gels  (10) . 

To  further  investigate  the  temperature  dependence  of  the 
substrate  on  the  metallization,  the  reservoir  was  cooled  by 
an  acetone-dry  ice  bath  to  180  K.  This  temperature  is  above 
that  of  the  glass  transition  of  the  polymers.  Samples  for 
System  B  were  metallized  at  this  temperature  and  the  resulting 
conductances  for  both  aluminum  and  indium  are  listed  in  Table 
2.3-VIII. 

The  values  of  the  conductances  were  intermediate  between 
those  at  300  K  and  those  below  the  glass  transition  temperature 
(T  )  for  both  metals.  The  results  for  indium  were  more 

o 

sensitive  to  the  substrate  temperature.  At  a  substrate  tempera¬ 
ture  of  300  K  the  indium  films  did  not  conduct,  whereas, for 
both  low  temperatures,  the  films  were  conducting.  There  appear 
to  be  several  competing  mechanisms  responsible  for  the  values 
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Table  2.3-VIII 

Conductances  of  A1  and  In  Films  Prepared  at  180  K 

System  Cure  Time(min.)  Conductances  (10~^  S  )  (a) 

_ A1 _ In  _ 

B  5  5  4 

10  7  14 

15  7  17 

a)  The  values  of  the  order  of  magnitude  are  reliable,  although 
there  was  spread  in  data  values. 
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of  the  conductances  which  include  the  primary  and  secondary 
nucleation  kinetics  (11)  and  metal  diffusion  in  the  crosslinked 
networks  (12) .  The  substrate  temperature  sensitivity  of  the 
conductance  indicates  that  a  nucleation  mechanism  is  predominant 
rather  than  the  glassy  state  of  the  substrate.  However,  further 
analysis  of  the  substrate  temperature  dependence  on  the  metallic 
penetration  into  the  polymer  is  needed  as  this  may  also  be  an 
important  factor  in  the  conductance. 

2.3.9  Conclusions 

The  inference  that  metallization  onto  a  glassy  rather  than 
a  rubbery  polymer  would  yield  a  metal  surface  with  improved 
electrical  and  optical  properties  is  supported  by  the  experimental 
results.  A  substrate  cooled  above  the  glass  transistion  temperature 
also  yields  better  conducting  films  than  those  of  an  ambient  sub¬ 
strate.  This  indicates  a  nucleation  mechanism  may  be  responsible 
for  the  substrate  temperature  dependence  of  the  film.  The  average 
reflectance  and  conductance  values  for  all  the  gels  prepared  by 
the  low  temperature  processing  were  higher  than  those  of  ambient 
processing  conditions.  An  improved  reflectance  is  also  found  for 
metal  films  with  low  temperature  processing.  This  is  in  agreement 
with  larger  and  smoother  sections  occuring  on  the  metal  surface  of 
these  gels. 

The  reflectance  and  conductance  values  for  the  metal  film 
processed  at  low  temperature  also  varied  with  the  metal  used, 
aluminum  or  indium,  the  silicone  system  and  the  curing  time. 

Among  the  gels,  the  aluminum  films  on  the  ultraviolet  cures  of 
System  B  yielded  the  largest  values  for  both  reflectance  and 
conductance. 
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2.4  FTIR  Studies  of  PDMS 


In  this  section,  the  IR  spectra  of  the  silicone  gel S  ystems 
A,  B,  and  the  Dielectric  Gel  taken  by  FTIR  are  presented.  The 
windows  of  the  first  set  of  spectra  were  MgF2  (supplied  by 
NVEOL)  and  had  a  cut-off  at  1000  cm--1-.  A  second  set  of  spectra 
were  taken  with  KBr  windows  to  allow  for  lower  wave  numbers. 

The  absorption  bands  were  correllated  to  the  silicone  polymer 
bond  structure.  This  enabled  the  kinetics  of  crosslinking  to 
be  followed  by  the  change  in  absorbance  of  reactive  bands.  In 
addition,  absorbance  measurements  were  taken  for  System  A  which 
yield  information  on  polymer  interactions. 

2.4.1  The  Spectra  of  System  A.  System  B.  and  Dielectric  Gel 

Samples  of  System  A,  System  B,  and  the  Dielectric  Gel  were 
coated  on  the  MgF2  windows  supplied  by  NVEOL.  The  systems  were 
crosslinked  for  the  following  exposure  times:  System  A  -  5  mins., 
System  B  -  10  mins.,  by  UV  ;  and  Dielectric  Gel  -  30  mins,  at 
413  K.  The  windows  were  found  to  transmit  between  5000  to  1000  cm-^ 
wave  numbers.  FTIR  spectra  were  taken  for  the  polymer  coated 
windows  over  this  region.  The  measured  transmission  spectra 
are  shown  in  Figs.  2.4-1,  2.4-2,  2.4-3,  for  Systems  A,  B,  and 
the  Dielectric  Gel  respectively.  Although  structure  in  the 
transmission  spectra  occurs  throughout  the  interval,  a  sharp 
absorption  is  noted  at  about  3000  wave  numbers .  This  is  a 
sharp  absorption  for  the  Dielectric  Gel  and  broadens  for  Systems 
A  and  B.  The  Dielectric  Gel  had  the  greatest  transmission  range 
as  compared  to  Systems  A  and  B. 

The  structure  of  the  spectra  and  relatively  low  transmission  for 
Systems  A  and  B  indicate  that  the  ultraviolet  active  side  groups 
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Infrared  Transmission  Spectra  for  Crosslinked  System 
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are  absorbing  in  the  IR.  If  this  is  so,  the  transmission 
spectrum  may  be  used  to  characterize  the  crosslink  density  and  , 
hence, the  modulus  of  the  microstructure.  In  order  to  separate  out 
the  silicone  bands,  we  have  identified  each  band  occurring  in 
PDMS  and  listed  them  in  Table  2.4-1. 


Table  2.4-1 

Bands  for  Determination  of  Crosslink  Density  ^ 


Bond 

Wave  Number (cm- 

Description  of  Band  or  Mode 

l 

Carbon 

C  -  C 

Methyl  Group 

1680  -  1620 

1750  -  1760 

Stretching;  Variable  intensity 

Stretching;  medium  intensity 

C  -  H 

2959 

C  -  H  stretch  in  SiCK-j 

Si  -  c 

800 

Si  -  C  bending  in  Si  (CH^^ 

ch3 

1377 

CH3  bending 

c  -  h2 

1263 

C  -  H2  rocking 

C  -  H 

1421 

C  -  H  bending 

Silicone 


Si  -  0 


1080  -  1020 


Intense  band  in  Polymers  of 
Silicone  and  Oxygen  due  to 
vibrations . 
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The  strong  absorption  for  the  silicones  at  3000  cm"^  can  be 
correllated  to  the  C-H  stretch  as  listed  above.  The  remaining 
modes  are  near  the  low  wave  number  cut-off  region  of  the  windows. 
In  Fig.  2.4-4,  the  transmission  spectra  for  dimethyl  and  phenyl 
substituted  silicones  are  shown .^There  is  very  little  structure 
in  the  region  of  5000  to  1400  wave  numbers  as  compared  to  our  TR 
spectra.  The  structure  is  believed  to  be  caused  by  crosslink- 
able  active  groups. 


2.4.2  Crosslinking  Kinetics  of  System  A 


By  extending  the  IR  measurements  into  the  low  wave  number 
region  with  KBr  window^  identification  of  the  crosslink  active 
side  group  bands  and  quantitative  measurements  of  crosslinks 
may  be  made.  We  have  therefore  retaken  the  IR  spectrum  of 
System  A  with  KBr  windows .which  allows  bands  in  the  lower  wave 
number  region,  from  1200  to  400  wave  numbers  be  observed.  We 
have  identified  those  bands  which  are  affected  by  ultraviolet 
cross linking. 

The  spectra  for  the  kinetics  of  crosslinking '  are  shown  in 
Figs,  2.4-5  to  2.4-7  for  which  a  System  A  silicone  gel  was 
consecutively  cured  for  1,  2,  and  4  minutes  respectively.  The 
scans  were  taken  from  4600  to  400  cm”^.  By  comparing  Figs, 

2.4-5  and  2.4-6,  we  notice  no  new  bands  occurring.  The  trans¬ 
mission  scale  for  Fig.  2.4-7  was  enlarged  for  ease  of  band 
determination  and  is  presented  in  Fig.  2.4-8.  We  notice  here 
that  at  about  1000  wave  numbers  there  is  a  sharpening  of  a  band 
that  does  not  occur  in  Figs.  2.4-5  and  2.4-6. 


Phenyl  tilicone  (film  on  sodium  chloride) 


Wavelength  in  microns 


Wavenumber*  in  cm’'1 


*  »  «  *  •  »  »  *  »•  II  il  l»  l«  IS 

Methyl-phenyl  silicone  (film  Winknfih  m  meroni 

on  sodium  chloride) 


Figure  2.4-4  Spectra  of  Typical  Silicone  Polymers 
From  Reference  1. 
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WAVENUMBERS 

Figure  2.4-5  Transmittance  of  System  A  cured  for  1  Min. 

The  six  bands  used  for  the  crosslinking  kinetics  of 
System  A  are  labelled. 


SYSTEM  A 
*  MINUTES 
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Fcr  the  1  and  2  minute  cures  the  transmission  mono- 
tonieally  increased  above  3000  wave  numbers,  whereas,  for  the 
4  minute  cure  there  was  a  monotonic  decrease  in  transmission 
for  this  region.  We  therefore  find  that  U.V,  curing  appears 
to  increase  the  opacity  of  the  gels  for  the  higher  wave  numbers 
on  the  border  of  the  near  infrared  (4000cm~^). 

In  order  to  determine  changes  in  the  bands  due  to  cross- 
linking,  six  reference  bands  were  chosen  as  marked  in  Figs. 
2.4-5  -  2.4-7.  According  to  Beer's  Law,  the  absorbance,  as 
defined  below,  is  directly  proportional  to  a  selected  bond  or 
concentration  as  follows: 

A  =  log  (1/1)  =  abc  (2.4-1) 

where  A  =  absorbance 

T  ■  transmission 
b  =  path  thickness 

a  =  frequency  dependent  cross  section 
c  -  concentration  of  substance  or  bond  of  interest 

By  forming  the  ratio  of  the  absorbance  of  a  non-reactive 
bond  (constant  concentration)  such  as  Si-CH^  or  Si  -0  to  a  U.V. 
reactive  double  bond  C  -  C,  the  kinetics  of  crosslinking  (open¬ 
ing  of  C  ■  C)  may  be  followed.  As  noted  from  Eq.  2.4-1,  this 
ratio  is  independent  of  sample  thickness. 

The  six  bands  identified  in  Fig.  2.4-5  are  listed  below  in 
Table  2.4-1. 
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Table  2.4-1  (Band  Identification) 
Peak  Number  Band 

1  SiO  band 

2  CH^  bending 

3  C  -  H  stretch 

4 

5 

6  C  =  C  band 


The  values  of  the  transmission  of  the  bands  as  a 

are  listed  in  Table  2.4-II. 

Table  2.4-II 

Transmission  Values  (T)  as  a  Function 

Feak(7„  Transmission)  12  3  4 

function  of  time 

of  Time 

5  6 

Time  (min) 

1 

19  23  27 

60 

64 

62 

2 

13  17  15 

58 

61 

40 

4 

64  68  59 

74 

67 

78 

The  absorbances  calculated  for  bands  1,  2,  3 

,  and 

6  are 

Table  2.4-III 

Absorbances  (-In  T) 

Band  No. 

Time (min) 

1 

2 

4 

1 

.721 

.886 

.193 

2 

.638 

.769 

.167 

3 

.568 

.823 

.229 

6 


376 


397 


107 


Forming  the  ratios  of  the  absorbance  of  band  6(C  -  C)  to  bands 
1,  2,  and  3, we  found  the  relative  change  of  the  carbon  double 
bond  due  to  crosslinking  as  listed  in  Table  2.4-IV. 


Table  2.4-IV 


Relative 

Absorbances 

of  Band  6 

Time (min) 

Relative  Band 

1 

2 

4 

3 

.662 

.482 

.467 

2 

.589 

.516 

.640 

1 

.521 

.448 

.554 

We  observed  that  thd  relative  absorbance  of  the  C  »  C  band 
to  the  C  -  H  stretch  band  (band  3)  decreased  for  the  2  minute 
cure  and  has  a  slightly  lower  value  for  that  at  4  minutes.  This 
appears  to  match  the  trend  of  the  previous  results  presented  for 
the  7.  extractables  for  System  A  in  the  bulk,  which  are  shown  in 
Fig.  2.3-1.  At  present  we  do  not  know  why  this  trend  does  not 
occur  for  bands  1  and  2. 

2.4.3  Crosslinking  Kinetics  of  System  B 

FTIR  spectra  were  also  taken  for  System  B  in  order  to  study 
the  kinetics  of  crosslinking.  The  sample  was  prepared  by  placing 
a  drop  of  System  B  on  a  KBr  window.  The  drop  was  successively 
exposed  to  U.V,  for  5,  10,  and  15  minutes.  The  values  of  the 
transmission  and  absorbance  spectra  were  measured  by  IR  for 
each  curing  time.  The  recorded  transmission  and  absorbance  spectra 
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in  order  of  increasing  cure  times  are  in  Figs.  2.4-9  to  2.4-16. 

A  list  of  the  bands  chosen  to  yield  ratios  of  relative  absorbance 
is  given  in  Table  2.4-V. 

The  absorbances  for  each  band  were  read  from  digital  output 
on  the  FTIR. 

Table  2.4-V 

Absorbances  for  Cures  of  System  B 
Absorbances 
Cure  Time  (min.) 


0 

5 

10 

15 

800.0 

2.00 

2.04 

1.88 

2.15 

1597.0 

0.08 

0.09 

0.10 

0.13 

1667.5 

0.07 

0.05 

0.00 

0.00 

1724.0 

0.00 

0.06  ’ 

0.11 

0.16 

2168.6 

1.03 

0.97 

0.83 

0.86 

2964.0 

1.37 

1.46 

1.27 

1.46 

2.4.4  Description  of  Spectra 

We  observe  from  the  Figs.  2.4-9  to  2.4-16  that  in  the  region 
of  1600  wave  numbers,  one  peak  at  1667.5  cm"1  decreases  with  curing 
time  and  another  at  1724  cm'1  increases.  The  carbon  double  bond 
stretching  occurs  at  about  1600  cm"1  for  conjugated  bonds  and 
as  high  as  1670  cm”1  for  the  nonconjugated  species.  The 
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Figure  2.4-13  Absorbance  of  Uncured  System 


Figure  2.4-14  Absorbance  of  System  B  Cured  for  5  Mins. 


2.4-15  Absorbance  of  System  B  Cured  for  10  Mins 
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Figure  2.4-16  Absorbance  of  System  B  Cured  for  15  Mins. 
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double  bond  transmittance  is  very  weak  in  the  unconjugated  case, 
but  very  intense  for  conjugated  cases  #  This  is  in  agreement 

with  the  value  of  1667.5  assigned  to  the  C  =  C  bond  as  the  PDMS 
structure  does  not  have  conjugation. 

We  also  note  that  for  wave  numbers  of  about  3350  there  is 
an  absorption  peak  developing.  This  could  be  due  to  water 
absorption  on  the  samples,  which  also  has  a  broad  absorption  at 
these  wave  numbers.  A  more  pronounced  effect  of  the  transmittance 
reduction  was  found  for  System  A  in  the  region  of  3000  -  4000 
wave  numbers . 

Ue  expect  that  the  bands  at  2964  cm“^  (C-H  stretch)  and  800  cm~l 
(9i-C  stretch)  remain  unchanged  during  the  crosslinking  reaction. 

In  order  to  trace  those  bands  active  in  the  crosslinking  kinetics, 
the  relative  absorbances  were  calculated  relative  to  the  C-H  and 
Si-C  bands.  This  method  of  ratios  normalizes  the  results  for 
different  thicknesses  of  the  polymer  film  scanned  by  the  IR. 

The  relative  absorbances  are  listed  in  Tables  2.4-VI  and  2.4-VII. 


78 


Table  2.4-VI 

Absorbances  Relative  to  C-H  Band 


Cure  Time (min.) 

Band 

(cm-1) 

800 

1597 

1667.5 

1724 

2168.62 

2964 

0 

1.45 

0.06 

0.05 

0.00 

0.75 

1.00 

5 

1.39 

0.06 

0.03 

0.04 

0.66 

1.00 

10 

1.48 

0.08 

0.00 

0.09 

0.66 

1.00 

15 

1.47 

0.09 

0.00 

0.11 

0.59 

1.00 

Table  2.4-VII 

Absorbances  Relative  to  Si-C  Band 


Sure  Time  (min.) 

Band  (cm"1) 

800 

1597 

1667.5 

1724  ' 

2168.6 

2964 

0 

0.04 

0.04 

0.04 

0.52 

0.80 

5 

1.00 

0.04 

0.02 

0.03 

0.47 

0.72 

10 

1.00 

0.05 

0.00 

0,06 

0.44 

0.67 

15 

1.00 

0.06 

0.00' 

0.07 

0.40 

0.68 

Except  for  the  data  for  the  5  minute  cure,  the  absorbances 
for  the  C-H  and  Si-C  bands  are  fairly  stable.  The  2168  band 
and  the  C  =  C  band  at  1667  decrease  during  the  curing  as 
expected.  However,  as  indicated  by  the  absorbance  data,  the 
crosslinking  kinetics  appears  to  saturate  after  5  minutes  of 
curing.  For  a  more  complete  study,  data  is  needed  within  this 
time  interval. 

2.4.5  The  Absorbance  of  System  A 

The  absorbance  in  the  infrared  as  a  function  of  the  silicone 
thickness  was  measured  for  System  A.  The  uniform  films  were 
prepared  by  using  mylar  strips  as  spacers.  The  thicknesses 
varied  between  6.4  pm  to  50.8  pm  and  are  listed  in  Table  2. 4- VIII. 

A  drop  of  silicone  was  deposited  on  a  KBr  window  between  the 
mylar  strips,  and  then  drawn  down  by  the  edge  of  a  microscope  slide 
guided  by  the  spacers.  This  drawdown  was  repeated  in  case 
polymer  relaxation  caused  the  thickness  to  increase.  The  films 
and  substrate  were  then  exposed  to  U.V.  for  5  mins.  This  is 
the  time  needed  to  obtain  near  maximum  curing  according  to 
sol-gel  content  analysis  presented  in  Section  2.3.3. 

The  IR  absorbance  curves  for  two  sample  thicknesses  of 
6.4  pm  and  23.4  pm  are  illustrated  in  Figs.  2.4-17  and  2.4-18. 
Although  the  absorbances  increase,  as  expected,  for  the  thicker 
sample,  there  appears  to  be  structure  in  the  peaks  for  the  23,4 
pm  sample.  This  indicates  that  new  bands  may  occur  when 
increasing  the  volume  of  the  gel.  These  jiew  bands  indicate 
a  polymer  -  polymer  interaction  in  the  gel  with 


Figure  2.4-18  Absorbance  of  23.4pm  thick  gel  of  System  A  Cured  for  5  Mins 
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increasing  volume.  If  this  occurs  we  may  not  expect  the  linearity 
between  absorbance  and  thickness  inherent  in  Beer's  Law  for  the 
data.  The  bands  for  this  analysis  are  those  that  are  not  altered 
during  the  crosslinking  reaction.  The  following  bands  should  be 
unreact ivej 


Wave  number 


(cm-1) 


Band  Assignment 


1022.00  Si-0 

1412.93  CH^  bending 

2964.70  C-H  vibration 


The  wave  numbers  and  absorbances  were  read  digitally  from 
the  FTIR.  The  data  are  listed  in  Table  2.4-VIII  below. 


Table  2.4-VIII 

Absorbance  of  5  minute  Cured  System  A 


Absorbance 


Sample  Thickness  (um) 

Wavelength  (cm 

1022.00  1412*93 

h 

2964.7 

6.4 

1.959 

0.280 

0,813 

12.7 

2.716 

0.398 

1.489 

19.1 

2.497 

0.363 

1.622 

23.4 

3.232 

0.637 

2.230 

36.1 

3.256 

0.860 

2.946 

50.8 

3.965 

0,915 

3.021 
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2. A. 6  Summary 

IR  spectra  over  the  interval  from  5000  to  400  cm"^  were  taken 

I 

for  System  A,  System  B,  and  the  Dielectric  Gel.  The  bands  active 
in  crosslinking  were  identified  according  to  the  bond  structure 
of  PDMS.  This  method  of  characterization  is  much  more  convenient 
;  than  that  of  swelling  for  the  case  of  thin  films.  For  long  exposure 

j  times  to  U.V. ,  the  silicones  are  altered  as  is  evidenced  by  the 

reduced  I.R.  transmission  at  the  longer  wave  numbers. 
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2.5  Deformation  Models 
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The  dynamics  of  thin  structures  have  been  modelled  using  a 
congressional  displacement.  The  first  model  which  used  an  electro¬ 
static  field  as  an  external  driver  was  solved  for  an  incompressible 
isotropic  medium  such  as  a  silicone  elastomer.  This  model  has 
been  extended  by  the  addition  of  the  more  correct  boundary  condi¬ 
tion  of  a  vanishing  displacement  field  of  the  structure  at  the  gel 
glass  substrate  interface.  This  boundary  condition  adds  a  geometri 
cal  shape  factor  into  the  dynamical  equations.  Equations  including 
the  viscoelastic  response  of  a  more  general  nature  are  presented 
in  Chapter  3.  The  equations  presented  below  have  been  used  to 
calculate  an  order  of  magnitude  of  static  deflections. 

2.5.1  Electrostatic  Deformation 

The  microstructure  can  be  envisaged  as  a  capacitor  with 
a  dielectric  of  silicone  elastomers.  A  diagram  of  the  structure 
and  the  labelling  coordinate  system  to  be  used  is  shown  in  Figure 
2.5-0.  The  elastomer  is  assumed  to  be  incompressible  and 
isotropic  in  the  unstrained  state.  The  deformation  due  to  the 
force  of  gravity  can  be  neglected. 

When  a  potential  is  applied  across  the  structure,  the  total 
energy  of  the  deformed  system  consists  of  two  terms,  the  electro¬ 
static  energy,  Eg,  and  strain  energy,  Eg ,  such  that 

F*ot)  =  Ee  +  Es  (2.5-1) 

For  a  uniaxial  compression  in  the  z-direction  the  net  normal 
stress,  Fz,  is 


F 


z 


?  E 

T7. 


£Ee 

T7~ 


3  E 

s 

T7~ 


(2.5-2) 
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where  the  partial  derivatives  are  evaluated  at  constant  voltage. 
At  equilibrium  the  net  normal  stress  is  zero  and  the  electrostatic 
attraction  is  balanced  by  the  mechanical  strain 

3  E  -  3  E 

_ e  a  _ s 

3Z  3  Z  (2.5-3) 

or  F  =  -F 
e  s 

For  a  parallel  plate  capacitor  with  dielectric  materials 
between  the  plates,  the  capacitance  is 


C  =  ee0A/d 

where  e  is  the  dielectric  constant, 

eQ  is  the  permittivity  constant, 

A  is  the  area, 

and  d  is  the  thickness  of  the  dielectric, 

The  stored  energy  Eg  in  the  capacitor  at  a  static  Voltage  is 
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Therefore  Fe  is 


F 

e 


EE 


(2.5-4) 


(2.5-5) 


From  the  statistical  theory  of  rubber  elasticity  the  strain  energy 
per  unit  volume  of  a  deformed  rubber  which  has  a  shear  modulus  G 
can  be  expressed  as 
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Es  °  7  G(X1  +  x2  +  X3  "  3)  (2.5-6) 

Where  Xp  x2,  X3  are  the  linear  extension  ratios  in  the  X,  Y,  and 
Z  directions  respectively.  For  a  incompressible,  isotropic 
uniaxial  compression  along  the  Z  axis,  where  is  the  initial 
length  and  d  the  deformed  length. 

^1^2X3  B  ^ 


Thus  .using  (2.5-6)  and  (2.5-7)  eq. 


F 


s 


G 


o 


(2.5-2)  becomes 


(2.5-7) 


(2.5-8) 


At  equilibrium, 


Thus,  from  (2.5-8)  and  (2.5-5), 


1  “oAv2 

'J~7~ 


G 


(2.5-9) 


For  unit  volume  AdQ  -  1,  and  the  deflection  for  a  potential  V 
becomes 


-  2  1 /3 
d  •  <do  -  do  “oV  /2R) 


(2.5-10) 


Equation  2.5-10  can  be  used  to  estimate  the  deformation  of  the 
micro8tructure  or  to  calculate  the  elastic  modulus  of  a  deformed 
gel. 
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It  was  shown  in  the  theoretical  derivation  (Equation  2.5-10) 
that  the  deformation  of  microelectronic  components  is  dominated 
by  three  factors:  the  applied  iroltage  ,  the  moduli  of  the  gels, 
and  the  microelectronic  component  geometry  (thiokness) .  The 
linear  percentage  deformation  ((d0  -  d)/dQ)  for  microelectronic 
components  under  various  static  voltages  has  been  calculated 
for  two  gels  with  shear  moduli  of  10^  Pa  and  10^  Pa.  In  the 
analysis,  5  pm,  10  pm,  and  20  pm  thick  samples  were  considered 
with  an  applied  voltage  between  50  and  400  volts.  The  results 
are  shown  in  Figures  2.5-1  and  2.5-2.  For  a  sample  with  a  thick¬ 
ness  of  10  pm,  450  volts  are  required  to  induce  a  107»  deformation 
in  a  gel  with  a  modulus  of  10^  Pa. 

Since  the  shear  modulus  G  (G  =  E/3)  for  System  B  obtained 
from  the  equilibrium  swelling  study  was  on  the  order  of  10^  Pa 
to  10^  Pa,  the  above  results  are  readily  applicable  to  estimate  the 
microcomponent  deformations. 

2.5.2  Deformation  by  a  Micro structure  with  One  Stationary  Electrode 

The  expression  derived  for  the  equilibrium  of  the  structure 
in  Section  2.5.1  did  not  implicitly  contain  the  boundary  condition 
ofja  stationary  or  fixed  surface.  For  the  micro structure  the 
bottom  electrode  is  mounted  on  a  rigid  glass  slide  and, therefore, 
at  the  interface,  the  gel  cannot  undergo  a  deformation.  A 
theoretical  treatment  for  small  compressions  including  this 
boundary  condition  has  been  derived  previously. 


APPLIED  VOLTAGE  (V) 


Figure  2.5-2  Component  Percentage  Deformation  Versus  Applied 

Voltage  (Shear  Modulus  x  lO^Pa)  From  Equation  (2.5-10) 
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The  overall  deformation  is  assumed  to  consist  of  the  super¬ 
position  of  two  simple  deformations: 

(1)  Homogeneous  compressions,  brought  about  by  the  same  deforma¬ 
tion  as  for  free  surfaces. 

Here,  the  force  is  given  as 


Fx  -  KEA  (dQ  -  d)  /d0 
where  E  «*  Young's  modulus 


(2.5-11) 


dQ  -  d  =  deformation  displacement 

K  «  a  constant,  which  is  equal  to  4/3  for  a  structure  where 
the  length  exceeds  five  times  the  width. 

(2)  A  second  force  due  to  the  outward  bulging  of  the  initially 
vertical  planes,  arising  from  the  requirements  that  the  volume 
of  rubber  is  conserved  and  that  a  stationary  deformation  exists 
along  one  surface.  This  force  is 

F2  -  (4/3)  S2EA(dQ  -  d)/dQ  (2.5-12) 

where  S  «  shape  factor,  depending  on  the  structure,  which  is 
equal  to  B/2dQ,  for  a  long  strip  structure. 

B  «  the  width  of  the  structure. 


The  total  force  Fg  is  given  by  the  sum. 


F  -  F,  +  F0 
s  1  2 

-  4/3  EA  (1  +  S2)(do  -  d)/dQ 
Using  the  equilibrium  condition, 


F«  "  “F» 
s  e 


(2.5-13) 


(2.5-14) 


end  Eq.  2.5-5,  we  find, 
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(1  +  S2)  EA  <dQ  -  d)/dQ  =  (3/8)  ee0AV2/d2 


(2.5-15) 


Figures  (2.5-3)  and  (2.5-5)  indicate  how  the  deformation 

depends  on  shape  factor  and  electric  field.  The  dependence 

of  the  displacement  on  applied  voltage  is  shown  for  two  values 

4  5 

of  the  modulus,  10  Pa  and  10  Pa.  The  curves  are  very  depen¬ 
dent  on  the  thickness,  dQ,  of  the  structure;  increasing  the 
thickness  reduces  the  electrostatic  force  and  .therefore, reduces 
the  deformation  whereas  decreasing  the  thifckness  increases 
the  deformation.  This  explains  the  trend  on  the  thickness  in 
the  figures  .  However,  increasing  the  thickness  also  diminishes 
the  shape  factor  which  hinders  the  deformation. 

When  the  displacement  versus  thickness  is  plotted  aB  in 
Figure  2.5-5,  it  is  found  that  there  is  an  optimal  thickness 
which  coincides  with  maximum  deformation.  This  optimal  thick¬ 
ness  depends  on  the  applied  voltage  and  the  elastomer  modulus. 


Figure  2.5-3  Component  Deformation  Versus  Applied  Voltage  (Shear  Modulus  x  lO^Pa) 
From  Equation  (2.5-15). 


't.OO  80.00  160.00  240.00  320.00  400.00  480.00  560.00  640.00 
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.gure  2.5-4  Component  Deformation  Versus  Applied  Voltage  (Shear  Modulus 
From  Equation  (2.5-15). 


d0  l/^m) 


Deformation  Displacement  Versus  Polymer  Thickness 
(Tensile  Modulus  x  lO^Pa)  at  several  Voltages 


97 


2 . 6  Experimental  Evaluation  of  Deformation 

An  optical  interference  experiment  was  set  up  to  detect  dis¬ 
placements  of  the  gels  to  a  fraction  of  the  wavelength  used  (HeNe 
laser  source) .  Problems  such  as  metal  breakdown,  charge  conductance 
through  the  gel  and  metallization  conditions  of  the  reflecting 
surface  were  investigated. 

In  order  to  insure  that  charge  was  stored  on  the  aluminum  elec¬ 
trodes  the  I/V  curves  were  measured  for  devices  of  several  electrode 
geometries.  Our  geometry  was  a  metal  film  for  the  bottom  electrode, 
the  other  electrode  was  a  metal  finger  pattern.  None  of  the  changes  in 
the  I/V  curves,  for  voltages  up  to  500  volts  indicated  breakdown. 

A  problem  arose  in  achieving  a  specular  reflecting  surface  for 
the  finger  electrodes.  This  led  to  investigation  of  the  optimuqi 
processing  conditions .of  the  metal  films  that  would  yield  good 
transmission  in  the  U.V.  and  still  be  conducting.  The  results 
indicated  the  optimal  choice  for  an  aluminum  film  to  be  120  to 
140  X  thick. 

2.6.1  The  Interference  Experiment 

A  schematic  of  the  system  is  illustrated  in  Fig.  2.6-1.  The 
laser  beam  is  divided  into  two  optical  paths  by  the  beam  splitter 
at  B.  One  path  uses  the  reflected  light  from  B  to  shine  on  the 
sample  at  C  then  through  B  to  the  screen  at  E.  The  other  path  is 
transmitted  through  the  splitter,  then  reflected  from  a  mirror 
at  D  and  combined  with  the  first  beam  at  E.  Care  was  taken  to 
eliminate  vibration  so  that , when  the  beams  overlapped  on  the  screen 
(E),an  interference  fringe  pattern  could  be  observed.  The  de- 
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formation  of  the  metal  surface  on  the  polymer  is  expected  to  shift 
this  pattern. 

2.6.2  The  Deformation  Experiment  Using  Gold  Electrodes 

It  was  found  that  a  sharper  interference  pattern  occurs  for 

a  gold  than  an  aluminum  film.  To  this  end,  a  structure  was  made 
of  a  glass  substrate  (0.846  cm  thick,  dielectric  constant  of  4.7), 
was  coated  with  a  silicone  dielectric  gel  (800  ym  thick,  dielectric 
constant  of  2,95),  and  then  a  thin  gold  film  of  400  &  was  evaporated 
onto  the  gel  at  cryogenic  conditions.  The  pattern  of  the  gold  film 
is  shown  in  Fig.  (2.6-2).  A  positive  potential  was  placed  at  the  T 
pattern  and  the  metal  holder  was  grounded.  A  separate  island  of 
gold  film  was  used  as  a  reference  which  did  not  deform  when 
placed  in  the  laser  beam  path  when  a  potential  was  applied  to  the 
T*  This  two-position  sampling  procedure  insured  that  the  shift  was 
caused  by  the  deformation  of  the  gel. 

Shifts  in  the  interference  pattern  were  observed  for  a  range 
of  voltages  applied  between  the  metal  holder  and  the  gold  film.  The 
deformation  of  the  gel  surface  was  calculated  from  the  measured 
shift  using  Eq.  (2.6-1) 

d  =  X  (AX)  / 2h  (2.6-1) 

where  d  is  the  deformation, 

x  is  the  laser  wavelength, 

AX  is  the  disdacement  of  the  fringe, 
and  h  is  the  width  between  two  constructive  interference  fringes. 


0  metal  holder 

®  GLASS  SUBSTRATE 
0  SILICONE  DIELECTRIC  GEL 
0  ft  0  GOLD  THIN  FILM 
®  LASER 
0  BEAM  SPLITTER 
0  MIRROR 
®  BEAM  EXPANDER 
©  SCREEN 


Figure  2.6-2  Schematic  of  the  Experimental  Set-up 
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The  values  of  the  deformation  for  h  =  2.4  cm  are  listed  in 
Table  2.6-1.  Since  the  potential  was  applied  across  both  the 
glass  and  the  gel,  the  correct  potential  across  the  gel  was 
calculated  using  the  thickness  dimensions  and  the  dielectric 
constants.  The  resulting  values  are  listed  in  the  second 
column  of  Table  2.6-1.  These  values  were  then  substituted 
into  Eq.  2.5-15,  the  theoretical  equation  for  the  displacement. 

The  deformation  is  plotted  for  different  values  of  the 
tensile  modulus  as  shown  in  Figure  2.6-3.  The  best  fitting 
value  of  the  modulus  is  about  3000  Pa.  This  is  a  low  value 
as  compared  to  the  value  expected  from  the  swelling  experi¬ 
ments  (10^  Pa.).  Experimental  errors  involving  the  measure¬ 
ment  of  the  thickness  of  the  gel  and  the  electric  field 
across  the  gel  may  account  for  this  discrepancy. 


2.6.3  Bulk  Current  Measurements  for  Aluminum  Electrodes 

Two  sets  of  capacitive  devices  were  constructed  for  the 
deformation  study.  One  set  consisted  of  an  aluminum  metal 
layer  on  glass,  a  crosslinked  20  ym  layer  of  silicone  (System 
B)  and  a  layer  of  indium  or  aluminum  metal  evaporated  when 
the  substrate  was  cooled  by  liquid  nitrogen.  The  second  set 
of  capacitors  had  metal  fingers  instead  of  a  metal  film  as 
the  electrode  on  the  glass.  The  fingers  were  20  um  wide 
with  a  separation  spacing  of  8  ym.  The  surface  for 
these  samples  was  indium.  In  order  for  the  deformation  to 
occur  and  be  maintained,  charge  must  be  stored  on  the  electrodes. 
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Table  2.6-1  -  Calculated  Displacements 


Applied  Voltage  Voltage  Across  the  Fringe  Deformation 

Polymer  Gel  Displacement  by  eq.(l) 

. V .  AX  Ccm>  d  C&) 

200  120  0.5  658 

250  150  .  0.7  924 

300  180  1:1  1455 

350  210  1.3  1709 

400  240  1.8  2373 

450  270  2.4  3164 

500  300  2.9  3822 


DEFORMATION 
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To  investigate  the  charge  storage  capacity  a  voltage  supply,  a 
capacitor  device,  an  electrometer  and  a  100  kn  safety  resistor 
were  connected  in  series  to  measure  the  I/V  characteristics  of 
the  devices. 

If  breakdown  of  the  gel  occurs,  the  current  versus  voltage 
curve  changes  slope  and  the  breakdown  voltage  can  be  deter¬ 
mined  where  this  change  occurs.  No  noticeable  current  was  measured 
in  the  tenth  of  a  milliampere  range  for  values  of  applied  voltages 
between  0  and  500  volts.  Since  breakdown  of  the  gel  did  not  occur, 
the  100  kn  safety  resistor  was  removed  and  currents  on  the  order 

_  Q 

of  10  amps  were  measured.  The  positions  of  the  leads  were 
fixed  in  place  since  minor  positional  adjustments  may  cause 
differences  in  readings  for  the  polymer  devices  for  such  small 
currents.  The  leads  at  the  voltage  supply  were  also  alternated 
so  that  the  I/V  curves  were  taken  for  both  polarities  of  the 
electrode.  The  I/V  curve  for  a  capacitor  is  shown  in  Fig.  2.6-4. 
The  open  circles  refer  to  the  bottom  electrode  on  the  glass  sur¬ 
face  being  connected  to  the  high  voltage,  the  dark  circles  for 
this  electrode  being  grounded.  The  I/V  curves  for  both  polarities 
appear  similar.  The  range  from  0  to  about  360  volts  appears 
linear  and  another  linear  section  of  increased  slope  lies  above 
360  V.  There  is  no  noticeable  saturation  for  either  polarity 
as  previously  reported  for  a  Ruticon  type  device. (1)  Figs. 

2.6-5  and  2,6-6  are  the  I/V  curves  for  the  second  set  of  structures 
which  have  metallized  fingers  for  the  bottom  electrode.  The  data 
for  such  a  structure,  Fig.  2.6-6,  should  be  read  with  care  as  two 
different  vertical  scales  are  used  depending  upon  the  polarity 
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Figure  2.6-5  I/V  Curve  for  an  Aluminum  Interdigital  Electrode 
on  System  B 
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Figure  2.6-6  I/V  Curve  for  an  Indium  Interdigital  Electrode. on 

System  B.  The  Vertical  Scale  on  the  Left  is  for  the 
Bottom  Electrode  at  a  Positive  Potential;  that  on 
the  right  is  for  the  Bottom  Electrode  at  a  Negative 
Potential. 
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of  the  voltage  applied  to  the  fingers.  In  both  cases,  when 
the  bottom  electrode  fingers  were  at  a  negative  voltage,  the 
current  was  largest.  The  I/V  curve  for  a  negative  polarity 
of  the  fingers  also  appears  to  have  some  curvature.  For 
positive  polarity  of  the  fingers,  the  I/V  curve  has  a  break 
in  slope  at  360  V,  with  a  sharper  slope  above  this  value. 

It  was  expected  that,  since  the  top  electrode  ha3  a  greater 
area  in  contact  with  the  polymer  than  the  fingers,  a  larger 
injected  current  would  result  upon  application  of  a  negative 
bias  to  this  surface.  However,  the  method  of  deposition  of 
this  metal  layer  is  different  than  that  of  the  fingers.  The 
effect  on  the  charge  injection  is  expected  to  be  different 
for  these  two  metallization  processes.  There  was  no  evidence 
for  a  breakdown  or  transport  of  a  large  quantity  of  charge 
for  either  sample. 

2.6,4  Surface  Metallization  Studies 

We  noticed  that  the  surface  metallization  of  the  capacitive 
devices  with  finger  electrodes  on  the  glass  substrate  was  not 
as  smooth  as  that  on  just  the  glass  substrate,  bight  could 
not  be  specularly  reflected  thus  causing  difficulty  in  form¬ 
ing  an  interference  pattern  to  observe  the  deformation.  The 
reflecting  metal  surface  on  the  gel  over  the  fingers  was 
studied  by  microscopy;  the  coarseness  of  the  surface  can 
readily  be  seen  in  Fig.  2.6-7. 

The  coarseness  of  the  metal  surface  was  thought  to  be 
caused  by  insufficient  or  non-homogeneous  crosslinking  thus 


creating  a  polymer  structure  with  surface  variation.  For  the 
formation  of  the  microstructure  element,  U.V.  light  must  pass 
through  a  thin  metal  film  and  then  enter  the  polymer  for  curing. 
This  film  must  also  be  conductive  to  act  as  an  electrode.  The 
deposition  process  was  investigated  to  optimize  the  thickness 
of  the  metal  film. 

A  major  characterization  parameter  of  the  film  is  its  thick¬ 
ness,  d,  which  can  be  estimated  from  equation  2.6-2  (2)  below: 

M 

d  «  - - e -  (2.6-2) 

*pt*  (1  +  <*/h)2)2 

where  Mg  -  the  evaporated  mass  (gms) , 

p  ■  the  density  of  the  film  material, 

h  -  the  distance  between  the  substrate  and  the  evaporated 
source, 

l  =  the  distance  from  the  center  of  the  substrate 


For  h  >>  i  and  an  Me  of  3mg.  of  aluminum  wire  a  6nm  metal  film 
layer  results. 

Previous  conductance  measurements  of  films  of  similar  thicknesses 
indicate  conductance  properties  that  are  not  reproducible.  The 
characteristics  of  the  metallization  of  the  glass  substrate 
is  not  only  controlled  by  the  thickness  but  also  by  the  kinetics 
of  film  growth.  Therefore  ,the  conductance  is  also  dependent 
upohthe  evaporation  rate  and  the  substrate  temperature.  An  experi_ 
mental  investigation  was  conducted  on  the  effect  of  evaporation 
rate  on  the  U.V.  transmission  and  conductances  of  the  films  in 
order  to  optimize  their  effectiveness. 


Ill 


2.6.5  Experimental  Details  for  Characterizing  Thin  Metal  Films 

Thin  films  of  high  purity  aluminum  were  prepared  by  thermal 
evaporation  from  a  tungsten  boat  in  a  high  vacuum  system  at  10"'* 
torr.  The  films  were  deposited  onto  glass  slides  at  room  tempera¬ 
ture.  The  rates  of  deposition  were  between  0.5  and  5.0nm/sec.  The 
film  thickness  was  measured  by  a  quartz  -  crystal  thickness 
monitor.  The  resistivities  of  these  films  were  measured  by 
a  four  point  probe.  The  U.V.  transmission  was  measured  by  a 
U.V.  visible  spectrophotometer. 

The  resistivity  of  thin  films  as  found  by  the  four  point 
probe  method  was  calculated  from  equation  (2.6-3)  (3) 

p/d  =  Rg  «  4.532  V/I 

(2.6-3) 

where  p  =  resistivity 

d  =  the  thickness  of  the  thin  film 
Rg  =  sheet  resistance 
V  =  voltage 
I  =  current 

The  calculated  resistivities  for  17  samples  are  listed  in 
Table  2.6-II.  For  some  of  the  samples,  a  variation  in  resistivity 
values  for  the  metal  surface  were  measured  and  for  these  cases 
the  term  "Random"  was  assigned.  The  values  in  Table  2.6-II  are 
also  plotted  as  resistivity  vs.  thickness  in  Fig  .  2.6-8. 

Table  2.6-II  shows  that  at  about  8nm  the  films  start  to  conduct, 
but  not  reproducibly  over  the  entire  surface. 
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Sample 

No. 

Thickness 

(A) 

Evaporation  Rate 

Resistivity 

(mr-cw) 

1 

46 

5-15  8/sec 

2 

55 

5-15  8/sec 

Random 

3 

59 

30  -  40  8/ sec 

7.75 

4 

66 

5-10  8/sec 

Random 

5 

88 

5-10  8/sec 

Random 

6 

109 

1-5  8/sec 

32,3 

7 

113 

10  -  20  8/sec 

1.54 

8 

129 

20  -  30  8/sec 

2.02 

9 

138 

10  -  15  8/sec 

1.98 

10 

149 

10  -  15  8/sec 

1.96 

11 

161 

10  -  15  8/sec 

1.13 

12 

165 

5-10  8/sec 

22.5 

13 

165 

10  -  15  8/sec 

3.17 

14 

181 

10  -  15  8/sec 

2.35 

15 

212 

10  -  15  8/sec 

1.13 

16 

215 

10  -  15  8/sec 

0.808 

17 

328 

10  -  15  8/sec 

0.838 

Comment 


1.  The" Random" 
means  there 
are  different 
resistivities 
at  the  different 
points  of  the 
sample. 


2.  Evaporation 
Rate  is  read 
from  thickness 
and  rate  monitor 


Table  2.6-II 

Resistivity  vs.  Thickness  for  Gold  Films 


Resistivity  (ufl 
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The  effect  of  evaporation  rate  on  the  resistivity  was  found 
by  comparing  the  values  for  samples  11,  12,  and  13.  The  film 
thicknesses  are  comparable  for  these  samples;  however,  the 
resistivity  of  these  films  is  inversely  proportional  to  the 
evaporation  rate.  This  result  is  in  agreement  with  thin  film 
theory  in  that  low  substrate  temperatures  and  high  deposition 
rates  enhance  the  formation  of  a  continuous  thin  film. 

The  U.V.  transmission  vs.  metal  thickness  appears  in 
Fig.  2.6-9.  The  cutoff  at  320  nm  is  due  to  the  properties 
of  the  glass;  this  restricts  the  usable  portion  of  the  U.V. 
spectra  to  between  300  to  400  nm  for  curing  the  polymer.  The 
U.V.  spectra  of  the  polymer  was  measured  to  determine  if  this 
range  overlaps  the  absorption  peaks  of  the  reactive  groups  for 
crosslinking.  If  the  frequencies  transmitted  by  the  glass  are 
causing  the  photchemistry ,  an  optimal  choice  of  an  aluminum 
metal  thickness  is  120  to  140  nm  deposited  at  a  fast  rate. 

The  resistivity  is  fairly  low,  1.0  -  2.0  yfl-cm,  and  the  U.V. 
transmittance  is  about  50%.  These  values  may  be  improved 
by  substitution  of  another  type  of  galss  that  transmits  over 
a  broader  U.V.  range  as  well  as  choosing  a  common  metal  (Au, 

Ag,  or  Cu)  that  has  a  higher  transmission.  Presently,  the 
above  investigations  are  being  extended  by  employing  a  special 
Nesa  glass  from  Corning  Glass  and  a  gold  metal  film  for  the 
surface  reflecting  electrode. 
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Chapter  3:  Microstructure  Dynamics 
3.0  Introduction 

The  dynamics  of  the  microstructures  depend  on  both  the 
elastic  and  viscous  properties  of  the  silicone  polymers.  Thus, 
equilibrium  compression  does  not  occur  instantaneously,  but 
depends  upon  the  frequency  of  deformation.  Similarly  the 
rise  time,  that  time  interval  which  the  structure  takes  to 
reset  itself  to  the  initial  isotropic  state,  does  not  occur 
instantaneously  but  is  dependent  on  the  viscous  para¬ 
meters.  At  kilohertz  frequencies  the  viscous  relaxation  cor¬ 
responds  to  localized  bond  rotations  of  the  polymers,  for 
lower  frequencies  in  the  Hertz  region, the  motion  corresponds 
to  whole  chain  motion  such  as  chain  slippage.  These  relaxations 
are  also  dependent  on  the  orientation  of  the  molecules  in 
the  gels  as  well  as  on  any  crystallinity  present.  The  above 
mentioned  effects  occur  for  the  bulk  of  a  structure  and  are 
intrinsic  properties  of  the  material.  However,  for  thin  films 
and  structures,  additional  complexity  of  the  relaxations  occur 
due  to  boundary  conditions.  The  boundary  conditions  confer 
relations  on  the  stresses  and  strains  that  couple  the  elastic 
and  viscous  loss  constants  and  may  also  involve  the  geometry 
of  the  sample.  Thus,  the  relaxation  times  for  modes  satisfying 
boundary  conditions  will  be  a  function  of  both  the  loss  and 
elastic  constants.  These  functions  for  plate  modes  have  been 
found  in  the  regions  of  low  frequency  and  long  wave  lengths 
and  are  presented  in  this  section. 

For  an  isotropic  material,  such  as  silicone,  there  are  four 
independent  elastic  and  viscous  constants.  The  plate  mode  re¬ 
lations  coupling  these  constants  may  be  used  to  determine 
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each  of  the  four  constants.  Once  these  fundamental  constants  are 
known  the  relaxation  times  of  any  geometrical  microstructure  with 
imposed  boundary  conditions  may  be  found. 

To  Illustrate  the  use  of  the  plate  mode  relations  in  deter¬ 
mining  elastic  constants , the  elastic  hysteresis  constants  for 
anisotropic  PVF2  were  evaluated. 


119 


In  general  oriented  polymers  films  are  anisotropic  and  the 
relations  between  the  elastic  and  loss  constants  are  found  ex¬ 
perimentally  from  acoustic  velocities  and  losses.  The  dispersion 
relations  for  the  losses  and  velocities  are  mathematically  complex 
for  an  anisotropic  plate.  However,  by  an  approximation  technique 
developed  by  Mindlin  (1),  the  dispersion  relations  for  the  low 
frequency  modes  may  be  directly  obtained.  This  procedure  was 
extended  to  viscoelastic  materials  and  the  equations  of  motion 
for  the  displacement  fields  were  obtained.  The  calculation  and 
resulting  solutions  of  these  equations  are  presented  in  Appendix 
3A.  In  the  following  section,  the  solutions  of  the  displacement 
fields  and  the  effects  of  the  boundary  on  the  coupling  of  these 
fields  are  illustrated. 


3.1  Description  of  Plate  Acoustic  Waves 


In  an  infinite  isotropic  medium  three  elastic  modes  of  pro¬ 
pagation  exist,  which  are  independent  of  each  other.  At  the 
free  boundary  of  a  plate  certain  of  these  modes  are  coupled.  The 
solutions  are  separated  into  four  classes  according  to  the 
directions  and  symmetry  of  their  displacements .  The  first  two 
classes  are  shear  waves,  i.e.,  one-dimensional  displacements 
of  layers  of  the  thickness  parallel  to  the  face  of  the  plate. 

The  two  shear  classes  depend  on  whether  the  symmetry  of  the 
displacement  is  odd  or  even  about  a  midplane  through  the 
thickness.  The  lowest  frequency  symmetric  shear  wave  is  illus¬ 
trated  in  Fig.  3.1-1,  along  with  the  lowest  frequency  waves  of 
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the  two  classes  to  be  discussed.  The  remaining  two  classes  of 
solutions  have  coupled  displacements  in  the  thickness  and  the 
propagation  directions.  The  lowest  frequency  wave  of  the  anti¬ 
symmetric  class  is  the  extensional  which  consists  of  antisymmetric 
displacements  in  the  thickness  direction  and  a  symmetric  longi¬ 
tudinal  displacement.  The  corresponding  lowest  frequency 
symmetric  mode  is  the  flexural  wave  whose  displacements  are 
symmetric  in  the  thickness  direction  and  antisymmetric  in  the 
longitudinal  direction.  These  modes  occur  in  the  region  of  low 
frequency  and  long  wavelength  and  are  of  experimental  interest 
in  characterizing  thin  films  as  to  their  elastic  and  loss 


constants . 
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3.1.1.  Tensor  Modulation 

Due  to  the  large  acoustic  losses  in  polymers  we  have  added 
to  the  elastic  constituitive  equations  which  yield  the  above 
three  (low  frequency)  modes  a  loss  term  which  represents  a 
viscous  or  a  hysteretic  loss.  The  dispersion  relations  derived 
for  the  above  three  modes  are  in  terms  of  the  components  of  the 
tensors  of  elastic  coefficients  as  well  as  those  of  the  viscous 
os  hysteretic  loss  teasor  .  As  a  general  example,  consider 
the  coefficients  compatible  with  the  symmetry  of  a  biaxially 
oriented  polymer  film  of  thickness  b. 

The  coordinates  are  choosen  such  that  and  are  in  the 
plane  of  the  film  and  x2  is  normal  to  the  plane.  For  this  film, 
the  orientations  of  the  amorphous  and  crystalline  regions  may 
not  align  along  the  draw  direction,  thus  a  triclinic  symmetry 
for  the  elastic  constants  is  assumed.  The  components  of  the 
elastic  tensor  conforming  to  this  symmetry  are 
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C14 

C15 

C16 

12 

C22 

C23 

C24 

C25 

C26 

13 

C23 

C33 

C34 

C35 

C36 

14 

C24 

C34 

C44 

C45 

C46 

15 

C25 

C35 

C45 

C55 

C56 

16 

C26 

C36 

C46 

C56 

C66 

(3.1-1) 


Subclasses  of  this  matrix  that  will  be  useful  are  those  of 


2mm  symmetry  for  PVF2  (13  constants)  of  isotropic  symmetry 
for  a  polymer  gel  (2  constants).  It  is  assumed  that  the  con- 


stituitive  constants  of  a  viscous  or  a  hysteretic  loss  will  con¬ 
form  to  the  same  symmetry  as  the  elastic  constants.  The  consti- 
tuitive  equations  for  the  two  loss  mechanisms  are: 


124 


pu. 


t  1 

C11 

I 

ul, 

11  + 

111  lit 

c15  u3 , 11  +  nll 

Ul.  11 

t  1 

n15 

11  " 

4^11- 

■  ^3 

1  1  t 

Cc51 

t  t 

c44 

c66 

b2 

T 

c16 

~ TT 

>  1 
(cll 

,'u2,UU  + 

lit- 

nll  u 

2,1111 

c66 

tit. 

+  nCi  u. 


(3.1-4) 


pu. 


"  C51  ul,ll  +  C55  u3 , 11  +  n51  ul,ll  +  n55  u3,ll 


t  t  » 

^54  c56  1,1  i  1 1  • 

T  jCl  rr]  Cc51  u2,lll»+  n51  u2 ,  llll3 
C44  c66 

»  t 

c56  ' *  '  *  •  i  • 

T  7^  Ccll  u2, 1111  +  nll  J2,llll] 

c66 


(3.1-5) 


For  the  flexure  mode  the  wave  equation  is 
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(3.1-6) 


The  constants  used  in  the  above  equations  are  defined  as  follows: 
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For  triclinic  symmetry  the  extensional  displacements  are 
coupled  to  all  three  displacements.  However,  a  decoupling  occurs 
for  the  2mm  symmetry  of  a  biaxially  oriented  polymer.  The 
elastic  constants  for  this  symmetry  are  given  by 
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(3.1-8) 


The  wave  equations  for  this  symmetry  decouple  and  become  for 
extensional  displacements 

Pul  “  cllul,ll  +  ni:ui,n 

Pu3  "  c55u3, 11  +  n55u3 , 11 


(3.1-9) 
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for  the  flexural  equation. 

When  the  function  for  a  wave  of  decaying  amplitude  is 
substituted  into  the  equation  of  motion  for  the  extensional 
modes,  the  real  and  imaginary  components  of  the  wave  vector 
are  given  as,  for  u^, 
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which  reduces  for  small  to 


'  2 

2  „//  nH  4 

a  =  P  /4  — r— *■  to 


The  characteristic  dependence  for  a  viscous  loss  of  a  on  w 
is  demonstrated  in  the  above  equation.  The  results  for  the  flexure 
mode  are  presented  in  Appendix  3A  and  will  not  be  discussed  in 
this  section.  The  dispersion  relations  for  a  viscous  loss  can  be 
converted  to  that  of  a  hysteretic  loss  by  the  substitution  of  the 

I  » 

hysteretic  coefficient  h^  for  n^w.  This  results  in  an  attenua¬ 
tion  coefficient  which  is  linearly  dependent  on  u  and  is  observed 
for  semicrystalline  polymers.  The  extensive  calculations  for 
the  hysteretic  loss  are  shown  in  the  Appendix. 

Further  reduction  occurs  for  the  evaluation  of  the  constants 
for  an  isotropic  medium  such  as  a  silicone  gel.  Here,  there  are 
two  independent  elastic  constants  and  C ^  and  two  viscous  con¬ 
stants.  The  coefficients  for  the  extensional  wave  equations  in 

terms  of  these  constants  are 
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for  the  displacement  Ug.  The  four  viscoelastic  constants  C^, 

^12’  nll’  an<*  ^12  can  ^*e  comPl-etely  determined  experimentally  from 
the  velocities  of  shear  and  extensional  modes  and  their  respective 
attenuations.  Experimental  results  to  illustrate  the  use  of  these 
equations  for  the  2mm  symmetry  of  PVF2  are  presented  in  the  next 
section. 


3 . 2  Experimental  Measurements 

The  measurements  were  made  with  delay  apparatus  following  that  of 
2 

Horvat.  It  consisted  of  two  blocks,  one  adjustable,  upon  which 
the  PZT  transducers  were  fixed  by  double- stick  tape.  The  delay 
line  is  shown  in  Fig.  3.1-2. 

The  PVF2  films  were  biaxially  oriented  Kureha  films  of  thick¬ 
nesses  25um  and  12ym.  Sample  strips  were  cut  1  cm  wide  along  the 
roll  direction  and  transverse  to  this  direction.  These  directions 
were  labelled  1  and  2  respectively.  One  side  of  each  strip  con¬ 
sisted  of  the  factory  finished  edge  which  was  aligned  along  the 
transducer  dimension  in  order  to  orient  the  PVF2.  No  further, 
analysis  of  the  film  orientation  was  made  and  it  was  assumed  that 
the  finished  edges  were  along  the  symmetry  axis .  The  films  were 
held  to  the  transducer  by  double  stick  tape.  Measurements  were 
taken  for  decreasing  values  of  the  distance  of  the  polymer  film  as  it  was 
found  that  removal  of  the  double  stick  tape  distorted  the  PVF2  film. 

For  the  propagation  of  the  SH0  model  (shear  mode) >the  PVF2 
strip  was  placed  along  the  length  direction  of  the  transducers 
such  that  the  displacement  motions  were  those  of  Ug.  The 
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exteneional  displacement  was  driven  by  aligning  the  strip  along 
the  width  direction  of  the  transducer.  (These  configurations  are 
illustrated  in  Fig.  3.1-3). 

Two  signals  were  generated,  one  a  pulse  mixed  with  a  cw, 
the  other  was  a  cw  mode.  As  shown  in  Fig.  3.1-4,  the  mixed-ipulse 
and  cw  were  fed  into  an  amplifier  before  being  sent  to  the  delay¬ 
line.  This  amplifier  was  necessary  since  6  to  12  volts  was  needed 
to  drive  the  transducers.  For  the  cw  measurement  the  signal  was 
directed  to  the  delay  line,  the  output  amplified  and  connected 
to  a  network  analyzer.  The  circuitry  is  shown  in  Fig.  3.1-5. 

3.2.1  Time  Measurements 

An  output  signal  in  response  to  the  input  is  shown  in  Fig. 
3.1-6.  The  input  pulse  width  was  kept  at  100  ps.  For  the  180 
KHz  signal  the  inout  voltage  was  6  volts  and  the  output  was  in 
the  millivolt  range.  The  delay  time  was  measured  from  the 
beginning  of  the  output  signal  to  the  beginning  of  that  of  the 
input . 

Measurements  were  taken  in  the  continuous  wave  mode.  The 
resonance  peak  was  located  on  the  network  analyzer  after  each 
sample  was  mounted  and  the  phase  of  the  propagating  wave  was 
measured  as  a  function  of  distance.  The  results  were  linear 
and  the  velocity  was  determined  from  the  slope.  Unfortunately,  in 
order  to  observe  the  phase  change,  the  velocities  could  only  be 
determined  with  a  precision  of  +  7%. 
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Figure  3.1-6  Input  and  output  for  an  SHq  wave  at  182  KHz 
for  25u  PVF2  film  (distance  =  1.9cm). 
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3.2.2  Attentuation  Measurements 

The  attenuation  measurements  were  made  by  measuring  the 
peak-to-peak  values  of  the  output  signal.  These  values  were 
normalized  to  the  input  voltage  and  the  decibel  loss  was  plot¬ 
ted  against  distance.  The  width  of  the  input  pulse  was  adjusted 
to  200  ys  to  insure  that  the  output  voltage  represented  the 
full  illumination  of  the  transducers.  An  attempt  was  also  made 
to  drive  the  input  transducers  with  a  continuous  wave  but  the 
output  voltage  did  not  monotonically  increase  with  decreasing 
distance,  because  of  the  reflection  of  the  wave  by  the  receiving 
transducer. 

3.2.3  Experimental  Results 

The  measurements  of  the  velocity  made  at  room  temperature 
on  PVF2  along  the  1  and  2  axis  are  listed  in  Table  3.1-1. 


Velocity  (m/s) 


Axis 

12pm 

25pm 

P„  Wave  1 

1728 

1897 

2 

1656 

1704 

Sh  Wave  1 

0 

1762 

841 

2 

1662 

893 

Table  3.1-1  -  Velocities  in  biaxially  oriented  PVF2  films.  The 


data  are  the  average  of  two  measurements . 
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Measurements  were  taken  on  two  thicknesses,  25  ym  and  12  ym. 
The  data  for  the  change  of  phase  angle  with  distance  for  the 
longitudinal  wave  along  the  1  axis  are  given  in  Figure  3.1-7. 
Resolution  of  the  velocities  could  be  made  between  the  1  and 
2  axis  for  the  25  ym  film  within  limits  of  the  experimental 
errors . 

The  attentuations  are  listed  in  Table  3.1-II  for  both 

shear  and  longitudinal  waves.  The  shear  attenuation  values 

are  in  agreeement  with  a  linear  frequency  dependence  to  data 

2 

taken  at  300  kHz  and  its  overtones. 


Attenuation 

in  dB/cm 

Axis 

12ym 

25ym 

PQ  Wave  1 

8.9 

10,6 

2 

4.4 

4.7 

ShQ  Wave  1 

9.7 

10.7 

2 

4.5 

6,2 

Table  3.1-II  -  Attenuation  at  200  kHz  in  biaxially  oriented  PVF2 
film.  The  data  are  the  average  of  two  measurements . 


The  plate  elastic  and  hysteretic  constants  according  to  the 
theory  derived  in  Equations  3.1-9  and  3.1-10  are  given  in  Table 
3.1-III  for  both  thicknesses  of  PVF2.  Scheyer  et  al  suggested 
that  the  hysteretic  loss  may  be  accorded  to  a  metastable  orientat- 
tion  of  polymer  chains. 
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Elastic  and  Hysteresis  Constants  in  N/m 


2 


12pm 

25pm 

C11 

5.34(109) 

6.44(109) 

C22 

4.91(109) 

5 . 19 (109) 

C66 (along  1) 

5 . 55 (109) 

1. 26(109) 

C66 (along  2) 

4. 94(109) 

1,42(109) 

"hi 

1. 50(109) 

2 . 37 (109) 

h22 

0.65(109) 

0.76(109) 

h66 (along  1) 

1 . 73 (109) 

0.21(109) 

h66 (along  2) 

0 . 67 (109) 

0. 14(109) 

Table  3.1-III  -  The  elastic  and  hysteresis  loss  coefficients 
as  calculated  by  Equations  17  and  18,  where 
~ —  _  „  Ci3  h33hi3Ci3 

U  11  '  ^  '  0^3 


and 


2 

Ci3h3i  +  h33Ci3 
C33  C33 


An  interesting  possibility  for  this  orientation  effect  is 
a  consequence  of  stress-induced  crystallization.  The  stress 
wave  could  stabilize  a  secondary  crystal  nucleation  site  on  the 
crystal  that  would  not  be  favored  by  thermodynamic  equilibrium 
under  ordinary  conditions.  Once  formed,  the  lifetime  of  such 
metastable  nuclei  may  be  enhanced  by  the  crystal  bulk  and  cause 
a  hysteresis  effect.  Thus  the  process  measurement  may  cause  a 
reorganization  of  the  sample  material. 
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Chapter  4:  Liquid  Crystal  Lens 
4.0  General 

Precision  optical  systems  are  limited  in  performance  by  a 
variety  of  phenomena.  Manufacturing  tolerances,  design  faults, 
material  instability,  and  thermal  expansion  all  contribute  to 
performance  degradation.  Optical  theorists  commonly  determine 
the  performance  degradation  in  terms  of  the  composite  effect 
on  the  optical  path  length  through  the  system,  especially  its 
variation  over  the  aperture.  Good  optical  performance  calls 
for  small  optical  path  length  over  the  aperture.  This  is  true 
not  only  for  variations  within  the  optics  but  also  for  variations 
in  the  external  propagation  path.  Thus,  when  viewing  a  star 
through  a  telescope,  the  star  image  which  should  be  a  small 
diffraction  pattern  may  be  spread  by  turbulence  produced 
variations  in  the  atmospheric  index  of  refraction.  If  these 
variations  could  be  compensated,  the  image  could  be  restored 
to  its  proper  size. 

Attempts  have  been  made  to  actively  compensate  for  the 
optical  path  length  variations  over  the  aperture  of  optical 
systems.  This  is  accomplished  by  servo  deforming  a  mirror 
in  the  optical  train  by  just  the  right  amount  to  correct 
the  path  length  errors.  To  achieve  a  significant  improve¬ 
ment  a  sizeable  number  of  points  on  the  mirror  must  be  con¬ 
trolled  leading  to  a  sizeable  number  of  servo  loops.  Further, 
the  servos  must  be  provided  with  a  command  signal  proportional 
to  the  desired  correction.  The  servos  must  have  a  bandwidth 
adequate  to  respond  to  the  temporal  variations  of  the  distur¬ 
bances  that  are  to  be  corrected.  In  the  field  of  adaptive 
optics,  systems  have  been  built  and  have  been  shown  to  per- 


140 


form  as  expected (1 , 2) .  A  significant  by-product  of  this  technique 
is  the  automatic  compensation  of  both  static  and  dynamic  structural 
path  length  errors  in  the  optics. 

Although  mirrors  have  been  the  primary  mechanism  for  adaptive 
optics,  any  mechanism  which  changed  the  optical  path  length 
(either  distance  or  index  of  refraction)  could  be  used. 

Lenses  play  a  critical  role  in  almost  all  optical  process¬ 
ing  systems.  From  the  crude  and  inexpensive  plastic  lens  to  the 
diffraction  limited,  ultra-expensive,  finely  polished  lens  made 
of  highly  pure  material,  their  role  in  the  optical  world  is 
unique.  Whatever  the  application  and  type  of  lens,  one  thing 
has  been  common,  the  need  for  focusing.  The  focusing  of  a 
lens  is  accomplished  by  mechanical  motion  of  one  portion  of  the 
lens  relative  to  the  desired  focal  plane  or  correspondingly 
to  other  elements  of  the  lens.  This  mechanical  movement,  even 
when  electronically  controlled,  requires  sophisticated  and  ex¬ 
pensive  hardware  to  be  precise. 

This  investigation  addresses  the  electronically  controlled 
creation  of  a  lens.  Since  the  actual  focusing  is  electronically 
operated,  it  can  readily  adapt  and  change  the  ultimate  focal 
point.  Also  since  the  index  of  refraction  is  controlled  point- 
to-point  on  the  lens ,  it  can  become  the  mechanism  for  adaptive 
optics. 

Unlike  electronic  focus  in  television  cameras  and  image 
intensifiers ,  it  is  not  the  bending  or  positioning  of  a  beam 
of  electrons  that  is  considered.  Instead,  it  is  the  direct 
electronic  control  of  the  index  of  refraction. 
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4.1  Background 

Several  methods  for  adaptively  changing  the  focal  length 
of  a  lens  have. been  investigated.  One  method  is  to  vary  the 
pressure  on  a  fluid  contained  in  a  lens  structure  C3) .  It  is 
limited  by  the  small  change  in  index  of  refraction  obtainable 
by  pressure  changes.  Another  technique  consists  of  putting 
rods  or  Kerr  cell  material  in  front  of  a  lens  (4) .  By  vary¬ 
ing  the  voltages  across  the  cylinders,  the  index  of  refraction 
and  ultimately  the  focal  length  of  the  entire  assembly  can  be 
altered.  Here,  the  limitations  are  granularity  of  the  cylinder 
stage  and  the  large  voltages  (kilovolts)  required. 

Another,  much  more  practical  technique  is  to  fill  a  hollow 
spherical  lens  structure  with  liquid  crystal  (5).  On  the  inside 
faces  of  the  structure,  transparent  electrodes  allow  varying  of 
the  voltage  across  the  cell  which  alters  the  index  of  refraction 
for  appropriately  polarized  light.  Limitations  here  arise  from 
constraints  on  cavity  size,  the  lens  must  be  small  for  uniformity 
within  the  liquid  crystal,  and  the  inapplicability  of  other 
than  a  single  incoming  polarization. 

Another  area  of  background  interest  is  the  small  lens  for¬ 
mulation  by  vortex  motion  within  a  liquid  crystal  cell  (6) , 
Although  not  of  practical  worth  as  a  lens,  it  does  demonstrate 
the  capability  of  varying  the  index  of  refraction  within  a 
liquid  crystal. 

4.2  The  Liquid  Crystal  Lens  Design 

Lenses  are  classically  created  by  varying  the  optical  path 
length  over  an  aperture  by  radially  shaping  a  medium  with  a 
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fixed  index  of  refraction.  As  Appendix  4A  shows,  if  the  structure 
has  parallel  walls  but  the  index  of  refraction  is  Varied,  the 
effect  is  the  sane. 

All  crystals  under  the  influence  of  an  electric  field  or 
applied  stress  exhibit  birefringence  (7).  That  is,  the  index 
of  refraction  can  be  varied  (See  Appendix  4B) ,  In  solids, 
however,  this  effect  is  small  since  only  distortions  of  the 
indicatrix  are  made.  Liquid  crystals  are  unique  since  the 
molecules  re-orient  under  the  influence  of  an  applied  field, 
either  electric  or  magnetic  (8) .  In  addition,  nematic  liquid 
crystals  are  uniaxial  and  the  optical  axis  is  coincident  with 
the  molecular  axis.  Thus,  as  the  molecules  re-orient  in  the 
presence  of  a  field,  the  entire  indicatrix  is  rotated  making 
large  changes  in  the  index  of  refraction  possible  (See  Appendix  4C) . 

Liquid  crystals  can  be  divided  into  two  groups.  If  the 
dielectric  tensor  is  such  that  the  component  along  the  molecular 
axis  (also  the  optical  axis)  is  greater  than  the  component 
perpendicular  to  the  axis  (e^  >  )  ,  the  crystal  is  said  to 

be  positive.  These  molecules  align  parallel  to  an  applied 
field.  If  the  reverse  is  true  (ejj  <  eJ,  the  crystal  is  said 
to  be  negative.  In  this  case,  the  molecules  will  align  per¬ 
pendicular  to  an  applied  field.  Either  type  can  be  used  for 
creation  of  a  lens.  Only  the  voltages  applied  to  the  electrodes 
will  differ  as  shown  in  Figure  4,2-1, 


To  obtain  total  polarization  capability  there  must  be  a  second 
stage  with  nematic  director  having  preferential  turn  in  Y- 
direction. 

Consider  the  case  of  a  positive  liquid  crystal.  In  order  to 
have  a  workable  lens,  the  molecules  would  be  in  a  homogenous  align¬ 
ment  without  any  voltages  applied.  Note  that  this  means  the  axes 
of  all  the  molecules  align  in  the  same  direction.  This  preferen¬ 
tial  direction  for  the  molecules  is  established  when  the  lens  is 
fabricated  and  can  be  accomplished  in  a  number  of  ways  (11,  19). 

The  simplest  method  is  preferential  rubbing  of  the  inside  surface 
prior  to  assembly.  As  the  voltage  is  applied  across  the  electrodes 
above  threshold,  the  molecules  tend  to  align  to  the  field.  At 
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saturation,  all  but  the  molecules  against  the  two  surfaces  will 
be  aligned  with  the  field  (See  Appendix  4D) .  If  the  preferential 
direction  is  in  the  x-direction  as  the  molecules  move,  the  index 
of  refraction  for  a  wave  traveling  in  the  z-direction  and  polarized 
in  x  has  changed.  The  index  of  refraction  for  a  wave  polarized 
in  the  y-direction  has  remained  constant.  Figure  4.2-2  shows 
graphically  the  alignment  state  and  the  corresponding  indicatrix 
location  of  n,  the  index  of  refraction. 

Thus  Figure  4.2-1  with  the  light  traveling  in  the  z  direction 
shows  a  configuration  which  with  appropriately  applied  voltages 
could  be  a  cross-section  of  a  lens  for  light  polarized  in  the 
x-direction.  To  account  for  a  y-polarized  component  a  second 
stage  would  be  required.  This  would  be  an  identical  stage  but 
rotated  90°  so  that  the  preferential  direction  is  now  in  the 
y-direction. 

It  has  been  assumed  in  the  previous  discussion  that  the 
electrodes  could  be  configured  to  give  the  required  radial 
variation  in  the  index  of  refraction.  During  this  analytical 
phase,  a  great  deal  of  time  was  spent  in  considering  the 
electrode  structure.  What  is  desired  is  a  radial  variation 
in  voltage.  It  would  appear  that  radially  symmetric  electrodes 
would  be  ideal  since  they  match  the  symmetry  required  of  the 
lens.  However,  they  require  difficult  mask  fabrication  and  any 
connections  made  to  them  create  obscurations. 

With  the  trouble  encountered  using  radially  symmetric 
electrodes,  rectangular  electrodes  similar  to  those  used  in 
liquid  crystal  displays  were  considered.  In  this  case,  it  is 
much  easier  to  address  individual  electrodes  and  mask  fabrication 
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Note: 


Figure 


4)  HOMOGENEOUS 


0  0  0 
0  0  0 
0  0  0 


b)  PARTIAL  ALIGNMENT 


0  0  0 
0  0  0 


C)  HOMEOTROPIC 


Since  angle  variation  is  in  X-Z  plane,  light  polarized 
in  Y-direction  will  be  influenced  by  nQ  in  all  three 
cases. 

4.2-2  Liquid  Crystal  Alignment  With  Light  Propagating 

in  the  Z-Direction  and  Polarized  in  the  X-Direction. 
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is  a  simple  process  but  the  device  inherently  does  not  match 
the  symmetry  of  a  lens. 

This  problem  was  solved  by  considering  a  square  aperture 
of  width  L  centered  at  the  origin.  Let  the  transmission  of 
the  aperture  be  given  by 


t(x,y)  =  rect  (£)  rect(^)  exp  r jkAn(x)  ] 


(4.2-3) 


Where 


k  =  wave  number 
A  =  thickness  of  the  material 

anc*  n(x)  ■  n  -  ■  =  index  of  refraction 

Just  past  the  aperture,  if  it  is  illuminated  by  a  unit-amplitude, 
normally  incident  plane  wave,  the  field  distribution  u^  is 


(x, y)  »  rect  (£)  rect  (~)  exp  (jka 


(n„  "  ) 


(4.2-4) 


Ignoring  the  constant  phase  term 


uL  (x,y)  =  rect  (£)  rect  (£)  exp  ("^-) 

But  note,  if  there  was  another  transmission  factor  directly  in 
front  of  the  aperture  with 

2 

n<y>  "  ne  '  hf 


(4.2-5) 


(4.2-6) 


Then 


t,(x,y)  tt2(x,y) 


“  rect  (£•)  rect  (£)  exp  ( 


-1k(x2  +  y2) 


(4.2-7) 


i 
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The  ideal  form  of  a  lens  is  obtained  with  linear  electrodes . 

By  cascading  two  stages  not  only  is  the  mask  fabrication  pro¬ 
cess  simplified,  the  required  symmetry  is  also  maintained. 

Figure  4.2-3  shows  a  single  3tage  cross-section  and  a 
representative  electrode  structure.  Figure  4.2-4  then  shows 
a  complete  lens .  "Complete"  means  that  it  works  for  all  in¬ 
coming  polarizations.  Note  that  four  liquid  crystal  stages 

A 

are  required.  Two  are  needed  for  x  polarized  light:  one 

A 

giving  the  required  variation  in  index  of  refraction  in  x  and 
one  giving  it  in  y.  Similarly,  for  the  y  polarized  component, 
two  stages  are  required. 
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Figure  4.2-4  Four  Stage  Liquid  Crystal  Lens 
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4 . 3  Liquid  Crystal  Lens  Performance  Analysis 

Although  the  lens  described  previously  possesses  the  radial 
symmetry  in  index  of  refraction  variation  of  a  common  circular 
lens,  it  does  represent  a  square  exit  or  entrance  aperture.  There 
is  no  loss  in  general  applicability  by  considering  the  performance 
of  the  square  configuration  since  either  a  round  or  rectangular 
field  stop  would  be  placed  in  the  lens  chain  whose  size  would 
become  the  limiting  factor  or  the  square  configuration  could  be 
the  desired  format  (for  example  in  television,  infrared,  or 
photographic  system  with  lxl  formats). 

The  common  measure  of  lens  fidelity  is  the  Modulation  Trans¬ 
fer  Function  (MTF) .  The  MTF  is  the  modulus  of  the  Optical 
Transfer  Function  (OTF) ,  which  specifies  the  complex  weighting 

factor  applied  by  the  system  to  the  frequency  component  at  (f  , f  ) 

x  y 

relative  to  the  weighting  factor  applied  to  the  zero  frequency 
component . 

To  evaluate  the  OTF  of  the  liquid  crystal  lens,  it  is 
first  necessary  to  examine  the  ideal  square  lens.  Ideal  here 
means  diffraction  limited.  An  imaging  system  is  said  to  be 
diffraction  limited  if  a  diverging  spherical  wave,  emanating 
from  a  point-source  object,  is  converted  by  the  system  into 
a  wave,  again  perfectly  spherical,  that  converges  toward  a 
point  in  the  image  plane. 

4.4  Diffraction  Limited  Square  Aperture 

Following  the  nomenclature  of  Goodman  (9) ,  the  OTF  of  a 
diffraction  limited  system  can  be  written  as 
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X  -  Wavelength  of  light 
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where 


and 


Substituting  equations  (4.3-3)  and  (4.4-4)  into  equation  (4.3-1) 
gives  the  OTF  as 


*<fx'fy>  = 


u  -  Ad.|fxi)a  -  Ad.|f j) 


lfx>  -  137 

|fyf  “  737 


(4.3-5) 


Let 


f  4  i. 

fc  I3J 


(4.3-6) 


f£  then  defines  the  spatial  frequency  along  the  axes  where  the 
OTF  goes  to  zero.  Using  the  definition  of  the  triangle  function, 
(9)  the  OTF  of  the  diffraction  limited,  square  aperture  can 
finally  be  written  as 


*cfx.y 


(4.3-7) 


4.4  OTF  of  the  Liquid  Crystal  Lens 

Due  to  the  electrode  structures  used  to  create  the  index  of 
refraction  variation  in  the  liquid  crystal,  a  diffraction  limited 
lens  will  not  be  formed.  The  wavefront  leaving  the  lens  will 
have  phase  distortions  caused  by  the  approximation  of  a  smoothly 
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varying  index  of  refraction  by  a  sampled  function.  Thus,  the 
liquid  crystal  lens  will  inherently  have  aberrations.  The  effect 
of  those  aberrations  on  lens  performance  as  measured  by  the 
predicted  OTF  is  the  intent  of  this  section. 

The  common  method  for  analytically  describing  aberrations  is 
to  model  them  as  a  phase-shifting  plate  within  the  aperture.  If 
the  phase  error  at  a  point  (x,y)  in  the  aperture  is  represented 

2ir 

by  kW(x,y)  where  k  ■  ~  and  W  is  an  effective  path-length  error, 
then  the  complex  transmittance  P  of  the  imaginary  phase-shifting 


plate  is  given  by 

F(x,y)  =  P(x,y)exp[ jkW(x.y) ] 


(4.3-8) 


where  P(x,y)  is  the  pupil  function  used  in  the  previous  section. 
The  complex  function  F  is  called  the  generalized  pupil  function. 

Again  following  the  nomenclature  of  Goodman  (9) ,  let  A(fxfy) 
be  the  area  of  overlap  of 

Xd.f  Xd.f  Xd.f  Xd.f 

P(5  +  — 2-*,  n  +  — and  P(« - n - j-^) 


When  aberrations  exist  the  OTF  can  then  be  written  as 


*d,f  xdif 

/  /  exp{jk[W(s+ — k — ,  n+ — *— 

a(fx,fy) _ ___ 

/  /  d?dn 
a(0, 0) 


.  f  Xd.f  xd.  f 

-w(e - n - i-^Jded 


(4.3-9) 


For  the  liquid  crystal  lens  described  previously,  since  the 
modulations  in  x  and  y  are  accomplished  in  successive  stages,  they 
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are  independent  and  equation  (4-9)  can  be  written 


Xd.f  Xd.f 

/  exp{jkCW(s  +  —i~£)  -  W(5 - H^)]}d? 

a(f  )  2  2 

*<f x-f v>  - - - T - 

X  y  /  d? 

a(0) 


Xd.f  Xd.f 

/  exp{ jk[W(n  +  — t-*)  -  w(n  -  -4-2)}d5 
a(f  )  2  2 

y 

- -  (4.3-10) 

/  dn 

a(0) 


Only  one  of  the  terms  (stages  in  the  lens)  needs  to  be  considered 
since  the  effect  in  the  x-direction  is  numerically  equivalent  to 
that  in  the  y-direction. 

Consider  the  electrode  configuration  depicted  in  Figure  4.4-5. 
The  electrodes  have  appropriate  voltages  applied  so  that  the 
index  of  refraction  is  identical  to  that  desired  at 


x  =  *(4n  +  l)a 
x  2 


n  “  0,1,2, 


N  +  1 
’  2 


and  N  is  the  number  of  electrodes.  Assuming  that  the  field 
underneath  each  electrode  is  uniform,  the  error  in  the  index  of 
refraction  under  the  electrode  is 


6n(x)  =  IDEAL  -  CREATED 


6n(x) 


Cnc  _  zfr3 


(4.3-11) 
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Figure  4.4-5  Liquid  Crystal  Electrode  Configuration  For 
Lens  Performance  Analysis 
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where  xQ  represents  the  center  of  the  electrode.  Figure 4.4-6  shows 
the  error  in  the  index  of  refraction. 

In  the  spaces  between  electrodes  it  is  assumed  that  cross¬ 
talk  exists;  i.e.,  the  spacings  are  small  enough  that  the  liquid 
crystal  molecules  are  reacting  to  a  field  strength,  the  magnitude 
of  which  varies  between  that  of  the  bordering  fields  established 
by  the  two  electrodes.  If  this  cross-talk  is  symmetric,  the 
index  of  refraction  would  also  be  appropriately  sampled  at 
x  =  ±  ;  n  -  1,  2,  .  .  .  The  dotted  segments  in 

Figure  4.6  illustrate  this  cross-talk  within  the  spacings. 

To  account  for  the  cyclic  nature  of  the  error  in  the  index  of 
refraction  and  the  increasing  error  with  distance  from  the  center, 
the  error  in  index  of  refraction  is  modelled  as 


fin(x)  =  -ar\ 


N2  -  (N  -  j)2 


]  cos  [2^] 
a 


(4.3-12) 


The  variation  given  in  equation  (4.3-11)  has  been  given  a  constant 
worst  case  value. 

Now 


(4.3-14) 


Using  this  expression  in  the  first  integral  of  equation  (4.3-10) 
yields , 

Xd.f  Ad. f 

N(f  )  -  /  exp{ jk[W(£+ — i— £)  -  W(5 - i-£)]}dS  = 

X  a(fx)  2  2 


1  Ad.f 

2  ir[5+— i-£] 

/  exp{-jka^[^][cos( - - — ± - ; 

a(f  )  2f  a 


nCC - i-S] 


.  cos  ( - - — ± - )  3>dc  (4.3-15) 


Using 


cos(a+g)  =  cosacosg  -  sinasing 


cos(a-g)  =  cosacosg  +  sinasing 


Equation  (4.3-15)  can  be  expanded  to 


9  N-'T’  _  7T  Xd  .  f 

N(f  )  *  /  exp{-jkazC-yr][^os^cos — — 

x  (f  )  a  4a 

x 


_  it  Ad.  f 

JLJSl 


sin  a sin  Za  ^ 


irAd.f  irAd.f 

(cos^cos — 2a'  X~  +  sin^sin — — ) ] } dg 


(4.3-16) 


Equivalently, 


CN-r] 


N(f  )  -  /  exp{jkaz  -r-  Csin^-i 


ir  Ad.  f 


a(fx) 


^sin_^]}dc 


(4.3-17) 
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For  simplicity,  let 

ka2[N-£]  lrXd.f 

A  -  - r_ sin  -  (4.3-18) 

So 

N(f  )  =  /  exp  {j  As  in  ^-}dc  (4.3-19) 

X  a<fx> 

This  expression  can  be  evaluated  by  expanding  the  integrand  in 
a  Taylor  series; 


2  3 

ejAsin0  3,  ^  +  jAsinQ  +  C  jAs^ine]  +  [ jAsinS  ]  + 


TT 


t r 


(4.3-20) 


and  integrating  term  by  term.  Considering  those  integrals  in 
order 


/  lde 
a(fx) 


d.  |f  | 
x 1  x1 

Z.2’.  _  .  de  “  i"  A dj^  |  f x | 


-A+  dilfx 


i 
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/  jAsin  2Li  d5 
a(f  )  a 


A  -  A  d  |  fx  | 

P 

-  *  +  Adi!fxj 

- 2 - 


jAsin  — -  dc 

CL 


-jA  f[cos  J  ( - V"^>  -  cos  7  (— 


t- Ad. if  I 
l 1  x 1 

T 


I  r  sln2^-  ds  -  /  d2si"2j^  d« 


a(f  ) 
x7 


-t+Xd.  f  . 
1  1  x' 

- 2 - 


(4.3-22) 


A2  a 


7  7 


-?  ~ 


?*-  ^dilfxl 

r-  sin^]  ^ 

4ir  a  .  ,  ,  ,  i  _  , 

£  -i  +  Xdi|fx| 

a  j - 


A2a,  it 


rcj«;ldil£xl)  -  ^(sinzt— 21-x'] 


Jt-Xd.  If 


-A+Xd. | f  | 

sin2[ - 

2 


_A2arTr 


«  w  A-Xd. lf  I 

r Ca(A_xdilfxl)  ■  sin2c - 2 — —  H 


(4.3-23) 


i 
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The  fourth  term  goes  to  zero  and  the  rest  of  the  series  is 
small  in  comparison.  Therefore,  the  numerator,  N(fx),  can  be 
found  by  substituting  (4 . 3-21) , (4 . 3-22) ,  and  (4.3-23)  into 
(4.3-19)  yielding 

.2  , 

N(fx)  -  U  -  AdjjfJ]  -  C&U  -Adi|fx|)  +  isin(|  £  -  Adjf^)] 

(4.3-24) 

Equivalently, 


A2a 


N(fx)  -  Cl  -  tj-Ki  -  Adi|fx|)  +  sin(|C£  -  Ad±  |  f x  |  )  (4.3-25) 


Since 


H(fx)  -  N(fx)//de  -  N(fx)/a 
a(0) 


(4.3-26) 


Then 


«(fx) 


(1  -  t^)- - Xf1-— ’ )  +  sin(J[A  -  Adi|fx|]) 


(4.3-27) 


Equivalently, 


A2 

H  fx)  -  (1  -  “ 


Adilfxl)i  a2  ^sinJU-AdJfJ]  (»-Adt|fx|) 
—————  +  -r—  C  _  "  3  » 


(4.3-28) 
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which,  once  again  using  the  definitions  of  the  triangle  and 
sine  functions,  can  be  written  as 


H( f  ) 


[1  - 


A2  .  A2 

X  +  T 


sine  (- 


^ilf> 


-)  ]  A 

c 


(4 


Figure  4,4-7  shows  a  comparison  of  two  configurations  of  liquid 
crystal  lenses  and  the  diffraction  limit  for  the  aperture  they 
fill. 


3-29) 
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Chapter  5s  Summary 

Our  basic  goal  has  been  to  evaluate  the  suitability  of 
polymer  systems  for  use  in  adaptive  optics.  Structures  such 
as  the  electrostatically  controlled  mirrors  ar.d  lenses  described 
in  this  report  require  the  mating  of  materials  relatively  new 
to  electronics  applications  with  microf abrication . technology . 
Microprocessor  technology  will  be  essential  in  electronically 
■calibrating"  structures  with  little  inherent  optical  quality. 

In  the  case  of  deformable  mirrors,  analysis  indicates 
that  elastomer  microstructures  should  be  capable  of  larger 
deformations  than  conventional  adaptive  optical  components. 
Microfabrication  creates  structures  with  large  area-to-volume 
ratios.  Thus  electrostatic  forces  can  overcome  inertial  forces. 

We  have  verified  that  such  materials  can  be  spread  in 
thin  sheets  and  formed  into  islands  smaller  than  20ym  on  a 
side.  The  techniques  developed  depend  on  using  elastomers  which 
can  serve  as  "resist".  Exposed  to  ultraviolet  light  through  a 
mask/  they  are  chemically  altered  and  the  unwanted  portions  can 
then  be  washed  away. 

We  have  also  demonstrated  that  conductive  and  reflective 
metal  electrodes  can  be  deposited  using  low  temperature  con¬ 
ditions.  The  thin  films  created  should  not  stiffen  the  polymer 
appreciably.  Indium  is  the  best  material  for  use  in  these 


structures . 
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Studies  of  the  elastomers  revealed  basic  physical  and 
chemical  properties  which  may  prove  useful. 

The  infrared  spectra  of  the  polymer  systems  were  correlated 
to  the  structural  bonds  of  the  silicone  polymers.  Using  this 
correlation,  the  kinetics  of  cross linking  were  then  followed 
by  comparing  the  absorbance  of  the  reactive  C-C  bond  to  the 
non  reactive  C-H  bond.  Only  for  polymer  system  (A)  do  kinetics 
as  found  by  IR  agreed  with  those  found  by  swelling  measurements. 

For  the  thicker  films,  the  absorbance  curves  indicates 
nonideality,  possibly  due  to  aggregation  of  molecules. 

In  addition  to  mirror  structures,  we  have  studied  the  con¬ 
cept  of  a  liquid  crystal  matrix  lens.  Appropriate  voltages 
will  produce  a  two-dimensional  variation  of  index  of  refraction, 
creating  a  spherical  phase  front  for  a  plane  wave  input.  Focal 
lengths  of  several  meters  appear  possible.  A  system  of  four 
tandem  nematic  lenses  has  been  shown  to  provide  focussing 
action  for  a  beam  of  light  with  arbitrary  transverse  polariza¬ 
tion.  With  a  dense  set  of  small  electrodes,  performance 
approaching  diffraction-limited  is  expected. 

Display  applications  cannot  tolerate  fringing  effects 
making  high  resolution  difficult  to  achieve.  For  the  lens 
application;  however,  the  fringing  effect  in  a  high  density 
array  should  create  a  graduation  in  index  of  refraction  greatly 
improving  phase  performance. 

During  the  next  year,  we  hope  to  make  continued  progress 
on  all  aspects  of  this  effort.  Experiments  are  planned  to 
expend  our  microfabrication  capabilities.  The  newly  developed 
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pattern  generator  will  be  used  to  write  directly  on  the  elastomer 
films  with  a  focussed  beam  of  ultraviolet  light  from  an  argon 
ion  laser.  Near  U.V.  wavelengths  will  be  obtained  by  special 
mirrors,  while  the  far  U.V.  will  be  introduced  by  doubling 
visible  light  in  non-linear  crystal.  This  technique  will 
free  us  from  the  need  to  subject  these  delicate  films  to  con¬ 
tact  "printing"  in  the  mask  aligner. 

At  higher  power  levels,  we  will  attempt  to  pyrolyze 
elastomers  which  do  not  crosslink  in  U.V.  light.  Thus,  new 
materials  may  become  available  for  microstructure  fabrication. 

Most  of  our  deflection  experiments  used  one-dimensional 
arrays,  rather  than  islands.  We  will  try  to  develop  techniques 
for  producing  bending  structures  to  produce  greater  deflection. 

We  expect  to  fabricate  a  lens  based  on  straight  nematic 
liquid  crystals  rather  than  the  commercially  available  twisted 
nematics.  Experimental  work  should  reveal  the  ultimate  limita¬ 
tions  to  be  oveicome  in  order  to  produce  a  true  image  forming 
lens . 
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APPENDIX  1A 

Static  Deflection  of  Microplates 


Let  us  estimate  the  requires  elastic  constant  of  the  elastomer 
and  the  magnitude  of  the  applied  voltage. 

The  torque  required  for  bending  the  elastomer  is  supplied 
by  the  forces  acting  between  the  charges  on  the  aluminum  films. 

The  maximum  torque,  fortunately,  occurs  at  the  foot  of  the 
elastomer  "L".  That  is,  it  occurs  where  the  elastomer  is  joined 
to  the  substrate  as  shown  in  Fig.  1.1-2. 

The  radius  of  curvature  p  at  the  foot  of  the  elastomer  "L" 
is  inversely  proportional  to  the  applied  bending  moment  M  and 
proportional  to  the  longitudinal  stiffness  constant  c  of  the 
elastomer  and  to  the  moment  of  inertia  I  of  the  elastomer  cross- 
section. 


Ic 

p  -  If 

The  deflection  angle  6  in  radians  of  the  elastomer  "L"  is 
approximately  equal  to 


(1A.1) 


e 


(1A.2) 


The  moment  of  inertia  I  of  a  rectangular  cross-section  with 
dimensions  b  and  w  about  an  axis  perpendicular  to  the  projection 
of  the  structure  shown  in  Fig.  1.1-2,  and  located  at  y  ■  0,  is 


I 


(1A.3) 
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where  w  is  the  dimension  perpendicular  to  the  projection  of  the 
structure  shown  in  Fig.  1.1-2.  Combining  equations  1A.1,  1A.2, 
and  1A.3  we  obtain  for  the  deflection  angle  e 


a  -  24M 

6 - T 

wb'u 


(1A.4) 


Let  us  next  estimate  the  bending  moment  M  due  to  the  force  be¬ 
tween  charges  on  the  aluminum  films,  assuming  the  two  aluminum 
films  form  a  parallel  plate  capacitor.  Granted,  this  is  a  fairly 
crude  assumption.  The  stored  energy  U  in  this  capacitor  at 
constant  potential  V  is 


(1A.5) 


The  capacitance  C  between  the  plates  is 


eoe  TT 


(1A.6) 


where  g  is  the  width  of  the  air  gap  and  e  is  the  relative  di¬ 
electric  constant  of  the  elastomer.  Here  eg  is  the  dielectric 
constant  of  vacuum.  The  dimensions  a  and  b  are  shown  in  Fig. 
1.1-2.  The  effective  force  F  on  the  upper  film  is  equal  to 
the  change  of  the  stored  energy  U  with  the  air  gap  g  at  constant 
potential  V 


F  -  -- 


(1A.7) 


Combining  equations  1A.5,  1A.6,  and  1A.7  we  obtain  for  the 


force 
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6enK  „2  „ 

0  V  &  r1_i_^n 

— ~  77  F  CI  +  F3 
b 


(1A.13) 


Assuming  a  deflection  angle  of  20°  (and  an  addition  deflection  of 
10°  at  the  knee  of  the  "L")  for  a  total  deflection  of  about  30°, 
with  a  value  of  5.106  volts  per  meter  for  V/b,  a  length  to  thick¬ 
ness  ratio  a/b  -  2  for  the  elastomer  arm,  and  a  relative  dielectric 
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constant  e  of  the  elastomer  of  about  4,  we  require  an  elastic 

/  2 

stiffness  constant  of  the  elastomer  of  about  1.5  .  1CF  N/m  , 
so  that  for  a  value  of  b  equal  to  5  urn  we  obtain  25  volts 
for  the  applied  voltage.  This  is  quite  reasonable. 

There  will  also  be  some  deflection  due  to  the  force  of 
gravity.  The  gravitational  bending  moment  is 

Mg  =agb2w[2b  +  ^  a]  (1A.14) 

where  a  is  the  density  of  the  elastomer  and  g  is  the  acceleration 
of  gravity. 

Calculation  yields  for  the  ratio  of  the  gravitational 

deflection  angle  0  to  the  electrostatic  deflection  angle  6 

6 


(1A.15) 


This  ratio  is  proportional  to  the  linear  dimension  of  the 
structure.  Thus,  the  gravitation  deflection  angle  would  be 
about  4°: 35'  for  a  device  with  b  =  lmn.  And,  of  course,  the 
required  voltage  is  also  proportional  to  the  linear  dimension. 
For  a  device  with  b  ■  1  mn  the  required  voltage  would  have 
to  be  5000V.  And,  of  course,  a  large  device  is  thus  more 
succeptable  to  acceleration  forces. 

One  further  note,  even  though  not  proven  here,  is  that 
the  selfresonant  frequency  is  inversely  proportional  to  the 
linear  dimension.  Also  the  response  time  is  linearly  propor¬ 
tional  to  the  linear  dimension.  Thus  a  5  vn  system  responds 
about  200  times  faster  than  a  lmn  system. 


APPENDIX  2A 

Experimental  Swelling  Procedure 
For  Crosslinked  Gel 


2A. 1  Materials : 

Spectral  grade  hexane 

Petri  dish  with  glass  plate  cover 

Planchette 

Balance 

Vacuum  oven 

2A.2  Crosslink  Density 

A  sample  of  the  network  should  be  weighed  dry  and  placed  in 
a  sufficient  amount  of  spectral  grade  hexane  such  that  it  is 
completely  immersed.  The  dish  is  covered  and  the  gel  is  swollen 
for  two  days;  it  is  then  removed  from  the  hexane  and  blotted 
twice  to  remove  excess  hexane.  The  gel  is  then  placed  in 
another  dish  and  cover  (preweighed)  to  determine  absorbed 
hexane.  The  temperature  should  remain  constant  at  25+l°C. 

2A.3  Sol  Content  (Uncrosslinked  Material) 

To  determine  the  amount  of  uncrosslinked  material  decant 
the  Supernatant  hexane  into  a  weighed  planchette.  Wash  the  gel 
and  swelling  flask  or  dish  with  hexane  again  and  pour  into 
planchette. 

Place  the  planchette  and  its  contents  into  a  vacuum  oven 
at  35°C  until  the  contents  are  evaporated  dry.  Reweigh  the 
planchette  and  determine  the  sol  fraction. 


APPENDIX  2B 

For mu 1 ations  of  Silicons  Gels 


The  materials  were  formulated 

in  the  following 

ratios  i 

System  A 

X2-7077 

97  pts.  by  weight 

X2-7078 

3  pts.  by  weight 

Sv3tem  B 

X2-7079 

75  pts.  by  weight 

X2-7080 

25  pts.  by  weight 

Dielectric 

Gel 

q3-6527  A 

50  pts.  by  weight 

Q.-6527  B 

50  pts,  by  weight 
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APPENDIX  2C 

Rheology  of  Spin-coating  Process 

The  rheology  of  the  spin-coating  process  has  been  studied 
by  Washo  (1)  based  on  the  assumptions  that: 

(i)  the  spinning  fluid  is  Newtonian; 

(ii)  no  evaporation  loss  of  the  solvent  occurs; 

(iii)  no  body  or  Surface  forces  affect  the  process. 

Washo  found  that  the  spinning  process  could  be  described  by 
the  relation: 

t  -  KCaFb  (2C-1) 

where  t  is  the  film  thickness,  F  is  the  rotational  speed  of  the 
disk,  C  is  the  volume  fraction  of  polymer  in  the  spinning  process, 
and  a,  b,  and  K  are  constants.  A  similar  investigation  of 
spin-coating  of  electronic  resists  has  been  done  by  Lai  (2) . 

The  log  of  the  thickness  of  microcomponents  (Table  (2C-I) 
versus  the  log  of  the  rotating  speed  F  is  plotted  in  Figure (2C-1) 
for  both  ultraviolet  curable  Systems  A  and  B.  The  average 
value  of  b  in  Equation  (2C-1)  obtained  from  these  samples  was 
different  from  the  theoretical  value  (-0.67).  This  difference 
may  be  due  to  the  fact  that  no  solvent  was  used  in  this  study, 
the  substrate  geometry,  and  the  viscoelastic  behavior  of  the 
polymers . 
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Table  2C-I 

Film  Thickness  VS.  Spinning  Speed 


Speed  (rpm) 


Thickness (yn) 

800 

1200 

1600 

2000 

2400 

2800 

3200 

3600 

System  A 

45 

32.5 

30 

25 

20 

17.5 

15 

12.5 

System  B 

40 

30 

20 

16.2 

13.7 

11.2 

10 

7.5 

THICKNESS  (/xm) 


speed 
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3A.1 


APPENDIX  3A 

Viscoelastic  Wave  Equations 


In  the  derivation  that  follows,  we  assume  the  constituitive 


relations  for  a  viscous  loss;  however,  except  for  the  boundary 
conditions  on  some  of  the  displacements,  the  form  of  the  dispersion 
relations  can.  readily  be  adopted  to  a  hysteretic  loss,  due  to 
the  symmetry  of  the  constituitive  constants. 

We  now  consider  the  generalization  of  the  equations  of 
motion  to  allow  for  a  dissipative  term.  The  equations  of  motion 
of  a  particle  are  given  by  the  extremum  principle 


6j  (T  +  W)dt  =  0 
1 


(3A-1) 


where  T  is  the  kinetic  energy  and  W  is  the  work. 

The  generalized  work,  W,  can  be  separated  into  two  terms 


w  -  EjTj5qj  +  EjQjSqj 


(3A-2) 


where  q^  is  a  generalized  coordinate,  T^  is  a  generalized  force 
that  is  derivable  from  a  potential  and  Qj  is  a  force  that  is  non¬ 
conservative  and  therefore  not  derivable  from  a  potential.  This 
force  is  the  source  of  the  loss  term  that  is  incorporated  into 
Hindlin's  theory.  After  performing  the  variations  with  respect 
to  generalized  modes  we  find  the  analogous  equations  of  motion 
for  a  continuum  to  be 

J  ^  (pUj  -  T^j  ^)6u,dVdt  “  J  ^  Qij  i<5UjdVdt  (3A-3) 

’  "i  ’  v  ci 
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where  p  is  the  density,  is  the  modal  displacement,  T^j  and  Q^j 
are  stresses  corresponding  to  conservative  and  nonconservative 
potentials  respectively.  These  equations  are  reduced  to  planar 
equations  by  assuming  the  three-dimensional  displacements  to  be 
written  in  a  power  series  in  the  thickness  coordinate  as 


uj  =  Enx2  uj (n)  (xl-  x3 ' 


(3A-4) 


The  superscript  on  the  displacement  refers  to  an  index  for 
uj •  By  judiciously  choosing  u..^n^,the  equations  of  motion  of 
the  three  lowest  frequency  modes  may  be  approximated  in  the  long 
wavelength,  low  frequency  limit.  These  modes  are,  in  increasing 
order  of  frequency,  the  flextrral,  the  straight-crested  thickness 
shear,  and  thickness  extensional  modes  as  shown  in  Fig.  3,1-1. 

It  i-s  fairly  obvious  that  a  combination  of  a  few  displacement 
fields  are  necessary  to  generate  these  solutions.  The  choice 
of  these  displacements  are  listed  by  order  index  and  are 
illustrated  in  Fig.  3A-1. 

The  zeroth  order  are  rigid  displacements  along  Xj ;  the  first 
order  displacements  contain  a  nodal  plane  and  are  antisymmetric 
with  respect  to  this  plane.  The  second  order  displacements  are 
symmetric  with  respect  to  the  nodal  plane. 


3A.2 


Stress  Equations  of  Motion 


Upon  substituting  the  power  series  into  the  variational 
equation  (3A-3)  f  and  integrating  over  the  thickness  the  remaining 
integrations  are  over  the  cross  sectional  area  and  time.  Since 
the  variations  are  arbitrary  the  integral  yields  the  equation 
of  motion 
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m  (n) 

Tij,i 


+  Fjn) 


p  B  u  • 
M  mn  j 


(m) 


-  n<!2jn'1)+  Pjn>) 

(3A-5) 


where  for  ■  1,3 


m  (n) 

Tij 


xn  T 

x2  ijdx^ 


b 

mJ  X2  Qijdx2 


(3A-6) 


for  i  =  2  integration  by  parts  yields 


rr.(n) 

T2j,2 


nT 


(n-1) 

2j 


Q2j . 2  =  Q  .(n)  -  nQ 


(n-1) 

2j 


(3a-7) 


where  and  »  are  the  nth  order  surface  traction  con¬ 

tributions 


where  F,(n)  x?  T 
J  -h  * 


and 


,  (n) 


2jdx, 


_^x2  Q2j 


dXr 


and  where 

rafn+1 

B  -  —2b - 

mn  (tn+n+1) 

-  0 


for  m  +  n  even 

for  m  +  n  odd  (3A-8) 


i 


l 
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Strain 

3A.3 

The  substitution  of  the  power  series  of  the  displacements  into 
the  strain  tensor 


Sij  "  %(ui,j  +  uj,i) 
yields 

-  %  tS<">  xj  <3A-9) 

where 


-(n) 

ij 


,<n) 


,(n) 


Ui» 3  +  Ui,i  + 


(n  +  1)  (6 


2jul 


(n  +  1) 


+  6 


2i  - 


(n  + 

j 


1) 


3A.4 


Constituitive  Relations 


The  conservative  stresses  are  given  by  Hooke's  law  as 


Tij  *  CijklSkl 


(3A-10) 


The  nonconservative  stresses  that  occur  in  the  variational 
principle  are  derived  in  the  case  of  linear  viscoelasticly  from 
the  Rayleigh  dissipation  function; 


T  "  %  nijklSijSkl 


■  1 


kl 


nijklSkl 


(3A-11) 


The  constituitive  equation  that  is  developed  by  combining 


both  stress  contributions  is 


184 


Cij  “  CijklSkl  +  n ijkl  Skl 


In  terms  of  the  planar  tensors  we  find 
/_  K  P  O  j_  _  C!  C^O 

cij  cijklbkl  n ijklbkl 


C3A-12) 


3A.5  Truncation  and  Reduction  to  the  Low  Frequency  of  Thin  Plates 


.  The  previous  formalism  can  be  used  to  the  first  order (n  =  1) 
and  boundary  conditions  can  be  imposed  to  reduce  the  set  of 
differential  equation  to  an  uncoupled  flexural  and  two  extensional 
modes . 


The  equations  of  motion  for  the  first  six  lowest  modes  are 
obtained  from  the  zeroth  and  first  order  approximation;  they  are 


•rCo)  +  r>^°)  +  _l  qCo)  =  2b  ou 

Tll|l  +  Qll,l  T31 , 3  q31,3  pu 


Co) 

1 


T  Co)  I  0Co)  i  -rCo)  .  pCo) 
U2.1  +  ^12,1  +  x32 , 3  +  -32,3 


2b  Pu<°> 


t>Co)  4-  f»Co)  +  tCo)  .  /-.Co)  _  »• 

T13, 1  +  Q13,l  +  *33,3  +  °-33,3  ”  2b  pu: 


Co) 


TC1)  .  0<1)  .  Td)  ,  0(1)  T0 

111,1  +  Qll,l  +  r31,3  +  g31,3  T 


21 


o  _  2bJ  ”C1) 

21  -  T  PU1 


TC1)  .  0C1)  ,  TC1)  .  0C1)  no  To  _  2bJ  Wl) 

T12,l  +  Q12,l  +  T32 , 3  +  Q32 , 3  Q22  T22  ~  ~7  pu2 

tcd  +  nd)  +  t(1>  +  n(1)  -  o°  t°  =  2b3,;;d) 

T13,l  +  qi3,l  +  T33 , 3  +  -33 , 3  q23  “  L23 


C3A-13) 
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The  displacements  are  governed  by  the  stress  resultants 

^and  ^  whereas  the  rotation  u£P  contain  contributions  from 
both  the  moments  T^P  ,  and  the  stress  resultants. 


The  Stress  Strain  Relations 

The  conditions  on  the  truncation  of  the  displacements  and 
boundary  conditions  for  the  first  order  equations  of  motion  are 


u<n>  -  0  n  >  1  i  -  1,3 

u|n>  -  0  n  >  2  (3A-14) 

The  tractions  u^=*  0  »  u^=«  u£2) vanish  at  the  surface  as 
well  as  the  following  stresses; 


t 


(1) 

22 


~  0 


At  the  first  order  the  kinet. -  and  potential  energies  con¬ 
tain  the  displacements  u^  and  u^^  .  The  equations  of  motion 
may  be  reduced  to  the  classical  extensional  equations  if  these 
displacements  are  not  coupled  into  them.  The  above  assignments 
of  vanishing  stresses  do  not  allow  these  displacements  to  appear 
in  the  energy  term  or  in  the  equations  of  motion.  Those 
boundary  conditions  result  in  equations  valid  at  low  frequencies 
since  and  u^  become  appreciable  for  thickness  modes  whose 

frequencies  are  greater  than  those  of  the  extensional  modes. 

The  zeroth  order  equations  become  upon  eliminating 
from  the  constraint  on 
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tll)  *  2b  tcUul?l  +  c13u3°3  +  c14  (“2°3  +  u3(1) 

+  C15  <u3?l  +  “$>  +  "u  +  ^ 

+  "ilul?l  +  "13U3?3  +  "14  +  u<2>) 

+  "i5  +  “$>  +  4  + 

‘M  *  2b  ^c31ul°l  +  c33u3°3  ♦  4  <u$  +  u<») 
+  C35  (u<°>  +  u<°>)  +  ^  <u<°>  +  u<b>) 

+  4u<b>  +  4^)  +  n;4<i(o)  +  -a)) 

+  n 35  <"#  +  4  -*  <u<°>  +  up>)3 

t23>  ‘  2b  Cc41ul°l  +  c43u3°3  +  c44  <u2°3  +  u3l>> 

+  c«  <u3?l  +  ul?3>  +  c46  <u2°l  +  “ll>> 

+  "41  “l?l  +  "«“3(?3  +  "44  <“2?3  +  “3l>> 

+  n43  (i<°>  +  i<°>)  +  n46  (i<°>  + 
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“  2b  Zc'5lu[°{  +  c53u3>3  +  C54  (u<°>  +  u^>) 

+  c55  (u3?l  +  ul?3>  +  C56  (u2?l  +  ul(1)) 

+  n51ul?l  +  n53u3?3  +  n54  (u2°3  + 

*12  "  2b  Cc61ul?l  +  C63U3?3  +  c64  (u 2 ?3  +  ^ 

+  c65  (u3°l  +  u$>  +  C66  (u2?l  +  Ul1>) 

+  n6l“l?l  +  n63“3?3  +  n  64  (i2?3  +  U3X>) 

+  n65  (u3?l  +  ul?3>  +  n66  (u2?l  +  ul1))] 


where  for  p,q  **  1,  2,  3 


Cp2C2q 


.  CP2n2g 


The  first  order  equations  are  obtained  by  eliminating  the 
three  strains  S^,  S^,  an(^  S^from  the  constraints  on  the 
stresses.  Upon  this  elimination  the  first  order  stress  displace¬ 
ment  equations  are 
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t(l)  _  2b3  rc,"ua)  +  c,,'ua)  +  c"’  (u<X>  +  U (!)) 
C11  “3”  CcU  U1  +  C13  u3  15  ^3,1  Ul,3; 


»fX>  +  <u'»Sl>  +  "15 ’  + 


«§>  -  £  + + <^n + 


«8>  ■  ¥  +  =;;'“3a)  ♦  •» 


+  4]l> 


(3A-16) 

Where  the  consecutive  eliminations  of  the  strains  yield  the 
constants : 


•  •  '  cp4c4c 

c  ■  c  -  *  i  ■  ■ 

Pq  P<1 


•  •  *  ' '  cp6c6q 

_ *  c _ _ trri 


c  m  c 

pq  pq 


and  viscoup  constants 


• "  1  cp4n  4q 

n  m  x\  -  ■  ■  ■  i  ■'* 

pq  pq  e' 


fit  1 


npq  "  npq 


(3A-17) 
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The  setting  of  t^  to  zer0  is  sufficient  for  u^^o  be  decoupled 
from  the  first  order  elastic  equations,  however,  it  is  observed 
that  the  vanishing  of  the  viscous  terms  coupled  to  u^is  accomplished 
by  setting  it  to  zero.  The  zeroth  order  stresses,  t^  and  t^, 
are  coupled  into  the  first  order  equations  of  motion.  These 
stresses  may  then  be  determined  in  terms  of  the  first  order  stresses 
and  inertia  terms  u^  and  u^  .  Another  coupling  of  zeroth  and 
first  order  displacements  arises  from  the  condition  that  the 
waves  of  propagation  be  straight-crested  along  the  thickness. 

This  imposes  the  restriction  that  S^°^and  Sg°^  are  due  to  dis¬ 
placement  components  of  a  rotation  and  vanish  through  constraints 


Yield 


-u<o) 
u2 ,3 


u(1)  -  -u(o) 
u3  u2 , 1 


(3A-18) 


When  the  above  restrictions  are  substituted  into  the  first 


order  equations  of  motion,  t^J^  anc*  *"23^  become 


t<0) .  [c;”u<»)n  +  3c;;'u<°>u  ♦  2^4%, 


(n  ci  +  2m  c)  u 


IS-'  u2,311  +  n53  u2,333 


•  (o) 

+  n 13  u2 , 331 


t  2n 


55  u2 , 331 


+  n 


11  u2 , 111 


] 


.(o) 

■23 


2b  Cc,”u(0)  +3c’,,u(0)  +  (c’"  +  2c',Su(o) 

Lc51  u2,111+  Jc53  u2 , 133  +  U31  +  C55  )u2,113 


c33  u2?333  +  n51  u2°lll  +  (n31  +  2n ** 


(o) 

55  '“2,113 


•  t  t  III.  ✓  \ 

(n53  +  2n35  )u2 1 1333 


(3A-19) 


The  rotary  inertia  terms  pu^^and  pu^^  are  assumed  small 
at  low  frequencies  and  not  included  in  the  above  expressions.  These 

two  equations  can  be  introduced  into  the  zeroth  order  stress- strain 
equations  by  solving  consecutively  for  and  from  the 

zeroth  ^rder  equations  in  terms  of  t^  and  t^°^  and  substituting 
into  the  remaining  equations.  The  zeroth  order  stress  displacement 
equations  are  now 


**!’  -  2b  +  c”'uc?>  +  c”’(u<”>  ♦  u<?>) 

+  n11  +  n13  +  n15  (u^°{  + 


+  c14  C1  '  c163  Ao)  ^  c16  .(o) 
“r~  +  c6 


'44  66 


'66 


t^-2b  +  +  + 
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1  1  t 

+  n31 

+  ^34 

C1  "  c363 

t  t 

t<°>  +  . 

c44 

c66 

t4  +-rr 

c66 

(o)  _ 


'31 


=  2b  [c 


51  ul,l  _r  ^53  u3 , 3 


(o) 


+  ct 


,(o) 


+  c 


55 


(u 


(o) 

3,1 


»$> 


(o) 


n  51  ui,i  +  n  53  u3 , 3 


(o) 


+  n 


55 


(u 


(o) 

3.1 


+ 


C54  [1  '  ^56]  t(o)  +  ^56 
c44  c66  4  c66 


(3A-20) 


The  above  three  stresses  are  those  remaining  after  the  re¬ 
duction  to  low  frequency  modes.  The  substitution  of  these  stresses 
into  the  zeroth  order  equations  of  motion  yields  the  develop¬ 
ment  of  the  displacements  due  to  general  triclinic  symmetry.  The 
resulting  displacement  equations  can  be  simplified  by  assuming 
an  infinite  dimensional  plate  in  the  direction.  The  displace¬ 
ment  equations  of  motion  for  this  case  are 


2bpUj 


(3A-21) 


Which  when  written  in  terms  of  the  displacements  become 
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Which  when  written  in  terms  of  the  displacements  become 


l  i  f 

b2  c14  rl  c  1 6 t  (o)  ,  "'-(o)  \ 

[1  _  i  (c51  u2>int+  n51  u2>11]j) 


b2  C11  ,  (o)  .  *  *'*(o)  v 

-  T~  <cll  u2,llll+  nll  u2,llll> 

c66 


•t°) 

PU3 


b2  C54  fl  C56-)rc,’'u (o)  +n",:(o)  , 

-  3~  -r-  Cl  -  -r-rJLc51  u2  mi+  n5l  ^2tlUr 

c44  c66 


b2  C56  rrM,,,W  +  n'"’<o>  n 

3 - "  tcll  u2 , 1111  +  nll  u2,im] 

c66 


bz 


u  rc  +  n ' ' 'i(o)  1 

3”  Lcll  u2, 1111  +  nll  u2, 1111 J 


(3A-22) 


The  first  two  equations  are  extensional  and  the  third  is  a 
flexural  equation.  The  viscoelastic  constants  couple  the  dis¬ 
placements  for  the  extensional  nodes.  The  effect  of  the  first 

2 

order  terms  are  observed  as  the  coefficients  of  b  .  It  is 
noted  that  the  flexure  mode  is  composed  entirely  of  first  order 
terms  and  that  it  is  uncoupled  from  the  extensional  displace¬ 
ments  . 


3A.6  Biaxial  Symmetry  Dispersion  Relation 


A  biaxially  oriented  film  of  PVF2  is  of  symmetry  2  mm  which 
has  elastic  constants  of  the  form 
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[Cij]  - 

C11 

c12 

c13 

0  0 

0 

C21 

c22 

c23 

0  0 

0 

C31 

c32 

c33 

0  0 

0 

0 

0 

0 

c44  0 

0 

0 

0 

0 

0  c55  0 

0 

0 

0 

0  0 

c66 

The  propagation 

equations  reduce  for  this 

symmetry 

pu<°> 

"  C11 

Jo) 

ul,ll 

1 

+  nu 

^(o) 

ul,ll 

pu<°) 

=  C55 

u(°) 
U3 , 11 

+  n55 

u(°) 

U3 , 11 

pu<°> 

b2 

~  '  1 

<CU 

u2°im+  m 

Jo)  \ 

u2,llll; 

(3A-23) 


(3A-24) 

When  the  function  for  a  decaying  wave  along  the  direction 
is  substituted  into  the  equation  of  motion  u^°\  the  real  and 
imaginary  components  of  the  wave  vector  are  given  as 


2  2 

k‘  -  pu4  Cl  +  (1  +  A^/c^  )%]/2cj1  (1  +  «V  /c^  ) 


11 


2 

a2  -  pa)2  [1  +  (1  +  u.2n21/c^1  )%3/2c^  (1  +  u,2n?1/c11  ) 


,  2 

11/C11 


which  reduces  for  small  to 


194 


2 

a 


P  /4 


(3A-25) 


The  same  dispersion  relations  occurs  for  the  shear  displace¬ 
ment  u<o)  and  are  found  by  substituting  the  constants  and 
into  the  above  equations. 

The  flexure  mode  yields  a  quartic  equation  in  a  and  k 
when  a  decaying  wave  is  substituted  in  the  equation  of  motion 
for  .  The  solutions  for  the  real  and  imaginary  parts  of 
the  wave  vector  are 


2 

a 


i  2  *2 

*|-(3p /c^b^)^  Cl  (1  -  4<d  n  ^  /cll  ^  )/2)  1 

2  2 

5  Op/c^jb2)*  a  -  ((1  +  (1  -  40)^^  /cu  )*/2)*] 


which  for  small  yields 

?  o  •  p  k  '  2  1  2 

a  ■  ti>  (3p  / c^b  )  n  ^CH 


(3A-26) 
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3A. 7  Hysteretic  Dispersion  Relations 

For  a  semicrystalline  polymer,  the  crystallites  act  as  cross¬ 
links  and  prevent  the  mechanical  motion  needed  for  a  viscous 
loss.  The  attentuation  has  been  evaluated  to  depend  linearly 
upon  the  frequency  which  is  indicative  of  a  hysteretic  loss. 

The  constituitive  equation  for  the  hysteretic  loss  is  of  the 
form 


c  s„  +  jh__s„ 
pq  q  pq  q 


(3A-27) 


where  h  is  the  loss  tensor  component  and  j  is-J-l.  It  is  noted 
that  there  is  a  direct  similarity  between  the  imaginary  component 
and  that  of  the  time  harmonic  viscous  loss.  An  analysis  following 
the  previous  derivation  leads  to  equations  of  motion  of  the  same 
form  as  the  viscous  equation;  however,  the  coupling  of  constants 
due  to  the  free  boundary  conditions  are  both  complex  and  are 

»  *  2 
C11  “  C11  "  c12  /c22  +  jh22 


.  * 

hll  "  hll 


c13h31^c22  +  ^h22 


(3A-28) 


The  extensional  equations  of  motion  for  u^0^  becomes 


CcH  -  (c122C22  +  h21h22c12^c22  +  h22  ^ ]  ul,ll 


+  j  Chjj  +  (ci22^22  ~  c12c22h21^c22  +  h22  ^  ul,ll 


-  pH 


(o) 

1.11 


(3A-29) 
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The  dispersion  relations  are  found  to  be 


2  pw' 


a(l  +  b2/ a2) 


HI  +  (1  +  b2/a2)^]/2 


2  _  p  < 


a(l  +  bVaz) 


C-l  +  (1  +  b2/a2)*]/2 


(3A-30) 


where 


a  "  C11  “  c12  c22  ~  h21h22c12 


C22  h22 


b  "  hll  +  c12  h22  "  c12c22h21 


c22  +  h22^ 


(3A-31) 


which  up  to  second  order  in  ^2,2  ^ecome 


a  *  c 


11  "  c: 


:12  c12h21h22 


b  "  hu  + 


c12  h22  c12h21 


(3A-32) 


Thus  the  dispersion  relations  for  the  hysteresis  loss  to  second 
order  terms  in  h^/c^j  are 


2,  •  2. 


(1  -  3/8  hu‘/cu‘) 


pu»2  C1  hll2 
2cll  ^  cll2 


(3A-33) 
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Similarly  the  dispersion  relations  for  the  remaining  modes 
can  be  found  by  substituting  into  the  viscoelastic  relations 
for  the  elastic  constant  and  for  the  viscous  constant 
For  the  viscous  case,  the  k  vector  decreases  due  to  second  order 
terms  in  the  viscosity.  Since  the  frequency  is  constant,  this 
yields  an  increase  in  the  phase  velocity.  A  similar  expansion 
for  the  hysteresis  loss  indicates  a  similar  trend  although  there 

is  also  a  second  order  decrease  in  the  elastic  constant  which 
tends  to  decrease  the  velocity.  This  decrease  is  a  second 
order  effect  and  does  not  appear  in  the  viscoelastic  treatment 
since  the  initial  condition  on  the  velocity  of  the  thickness 
displacement  u^  was  set  to  zero  for  the  viscoelastic  case. 

3A.  Summary 

The  displacement  equations  of  motion  of  the  three  lowest 
frequency  modes  of  propagation  in  anisotropic  lossy  plates  have 
been  derived  by  Mindlin's  formulation.  The  boundary  conditions 
of  vanishing  stresses  across  the  plate  couple  the  loss  constants 
to  the  elastic  constants  in  the  equations  for  the  remaining 
stress  equations.  The  dispersion  relations  derived  using  these 
stresses  yielded  the  attenuation  dependence  for  the  extensional 
waves  to  first  and  second  powers  of  the  frequency  for  the 
respective  hysteretic  and  viscous  losses.  However,  the  flexural 
attenuation  for  a  viscous  loss  was  to  the  three- halves  power 
of  the  frequency. 

The  coupling  of  the  equations  are  symmetry  dependent  and, 
for  a  biaxial  film  of  2  mm  symmetry  such  as  PVF2,  the  extensional 


shear  wave  viscoelastic  constants  are  uncoupled  and  can  be  evaluated 
directly  from  experiment.  The  elastic  and  hysteretic  plate  mode 
constants  along  each  principle  axis  were  evaluated  including  the 
values  of  the  uncoupled  shear  constants. 
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APPENDIX  4A 

Creation  of  a  Thin  Leris  By  Index  of 
Refraction  Variation 

Consider  the  geometry  shown  in  Figure  4A-1  where  a  material 
of  width  A  and  variable  index  of  refraction,  n(x,y),  is  centered 
on  the  z-axis  and  lies  in  the  x-y  plane.  A  wave  at  coordinates 
(x,y)  will  suffer  a  phase  delay  in  passing  through  the  material 
given  by 

0(x,y)  =  kAn(x.y)  (4A-1) 

Where 

k  »  —  =  Wave  Number 


X  =  Wavelength 


In  attempting  to  vary  the  index  of  refraction  to  match  the 
spherical  curvature  of  a  thin  lens,  assume 


n(x,y)  =  nQ  I  1  - 


>i  ~Vri  - 


2  2 
x-y 


)+h^K 


2  2 
x-y 


(4A-2) 


Where  Figure  4A-2applies  for  definitions.  The  phase  delay  then 
becomes 


0(x,y)  =*kA  nQ  1 


(4A-3) 


Equivalently, 

0(x,y)  -  kAnQ  -  kO^-DRj. 


(4A-4) 
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Using  the  Paraxial  Approximation 


x2  +  V2 


(4A-5) 


Equation  (4A-4)  becomes 


0(x,y)  «  kAn  -  k(n  -  1)R,  1  -  1  - 


or  finally, 

2,2 

0(x,y)  =  kAnQ  -  k(nQ  -  1)  X 


2  ,  2 
x  +  y 


(4A-6) 


(4A-7) 


The  Corresponding  Transmission  Function  is 


2  2"| 

-ft-(x.y)  -  exp  j0(x,y)  =  expjjkAnJ  exp  [-  jk(nD-l)  ^-J  (4A 

But  using  the  definition  of  focal  length  for  a  plano-convex  lens 


1  4  "O  -  1 

T  =  R, 


(4A-9) 


Equation  (4A-8)  becomes 


it(x,y)  -  expfjkAnJ  exp  (x2  +  y2)] 


(4A-10) 


Equation(4A-10)  is  the  transmission  function  of  a  thin  lens 
Note  that  n(x,y)  could  have  been  considered  to  increase  rapidly 
and  if  appropriately  chosen  would  have  resulted  in  a  plano-con¬ 
cave  equivalent  lens. 

The  form  of  the  index  of  refraction  given  in  equation  (4A-2) 

2  2 

can  be  simplified.  Letting  r  -  x  +  y  ,  it  becomes 


n(r)  -  n. 


<"o  -  l> 


i  -Vri  - 


(4A-11) 
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Equivalently, 
n(r)  -  no 


Again  using  the  paraxial  approximation 


(n  -  l)Ri 
n(r)  -  n^  -  -2__ - l 


1  -  1  +  ±-~ 
2Rf 


(tu  -  l)r' 
n(r)  -  n  -  0 


2aR, 


(4A-12) 


(4A-13) 


Finally,  using  the  definition  of  focal  length  given  in  Equation 
(4A-9)the  required  variation  of  the  index  of  refraction  becomes 


n(r) 


(4A-14) 


( 
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APPENDIX  4B 


Birefringence 


Crystals  generally  are  electrically  and  thus  optically  anisotropic.  An 
electric  field  produces  a  polarization  in  a  crystal  which  is  not  always 
a  scalar  times  the  field.  Instead  the  magnitude  and  direction  of  the 
polarization  depend  on  the  direction  of  the  applied  field  in  relation  to 
the  crystal  lattice.  As  will  be  shown,  one  important  consequence  of 
this  phenomenon  is  that  the  velocity  of  light  propagation  within  the 
crystal  is  a  function  of  the  direction  of  propagation  and  the  polariza¬ 
tion  of  the  light.  • 

A  monochormatic  ray  of  light  propagating  through  an  anisotropic  crystal 
is  split  into  two  rays.  These  orthogonally  polarized  rays  generally 
travel  at  different  speeds  and  in  different  directions.  This  phenomenon 
is  known  as  birefringence  or  double  refraction.  Natural  double  refraction 
is  exhibited  by  20  of  the  32  solid  crystal  classes  and  all  liquid  crystals. 


In  an  anisotropic  material,  the  relation  between  the  displacement,  D,  and 
the  electric  field  intensity,  E,  can  be  written 


Di  "  fEoeijEj  1  ~  1’2’3  <«-!) 

where  £  «*  Free  space  permittivity 

o 

£.,  *  Components  of  the  dielectric  constant  tensor 


C4B^2) 
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Where  Einstein  notation  -  summation  is  to  be  performed  over  that  index 
which  appears  gore  than  once  on  any  term  (in  this  case  j  »  1,3,2)-  is 
used.  Equation  4B-2  written  out  in  matrix  form  is 


D 


1 


1 

e 

o 


t 

e 

c 


11 

C12 

E13 

21 

E22 

E23 

31 

e32 

E33 

C4B-3) 


The  single  quantity  e  in  the  isotropic  case  is  replaced  by  an  array  of 
nine  components  which  relate  the  three  components  of  D  to  the  three 
components  of  E.  D  is  no  longer  necessarily  in  thq  same  direction 
as  E. 


The  Maxwell  equations  are  still  , 

V  X  E  -  ^§f  (4B-4) 

and  t 

V  x  H  *  J  +  3D  (4B-5) 

3t 


Transparent  crystals  have  low  conductivities  and  are  not  ferromagnetic.  Thus 


J  -  0 

and  B  -  PpH 

Using  these  equations  (4B-4)  and  (4B-5)  become 


V  x  E  - 

'VQ 

3D 

at 

and 

V  x  H  - 

(4B-6) 

(4B-7) 


(4B-8) 

(4B-9) 


! 
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0 


(4B-13) 
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2 

This  is  a  quadratic  equation  in  v  .  If  the  negative  values  of  v  are 
discarded,  there  are  in  general  two  values  of  v  associated  with  each 
wave  normal  1.  The  E  aud  D  fields  associated  with  each  value  of  v  are 
found  to  be  linearly  polarized.  Further,  D  vectors  with  different 
velocities  are  orthogonal  to  each  other. 

The  relative  vector  configuration  is  shown  in  Figure  4B-1.  In  an  iso¬ 
tropic  material  or  free  space,  D  is  collinear  with  E  and  the  Poynting 
vector,  ?  =  E  x  H,  is  collinear  with  T.  In  an  anisotropic  material 
D  is  related  to  E  through  a  tensor  relationship  and  is,  therefore, 
not  necessarily  parallel  to  E.  It  follows  that  S  is  not  necessarily 
parallel  to  T.  This  means  that  the  direction  of  energy  propagation  is 
not  the  direction  of  phase  advance. 

• 

The  value  of  c/v  for  each  of  the  two  possible  waves  is  defined  as  the 
refraction  index,  n  for  that  wave.  The  refractive  indices  of  the  two 
waves,  as  functions  of  the  direction  of  their  common  wave  normal, can 
be  obtained  by  drawing  the  indlcatrix,  the  ellipsoid  shown  in  Figure 

4B-*. 

If  Xj^,  x2,  and  x3  are  the  principal  axes  of  the  permittivity  tensor, 

the  indicatrix  is  defined  by  the  equation. 

2  2  2 

^1  +  X2  +  X3  -  1  (4B-14) 

2  2  2 

nl  n2  n3 


_ s _ 


ISOTROPIC  ANISOTROPIC 


Figure  4B~1 

VECTOR  COMPARISON  OF 
ISOTROPIC  AND  ANISOTROPIC 
MATERIALS 
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where  n^,  nj,  and  n^  are  the  principal  refractive  indices.  The  lndicatrix 

has  the  following  important  property.  Draw  through  the  origin  a  straight 
line  OP  in  an  arbitrary  direction.  Draw  the  central  section  of  the  indic- 
atrlx  perpendicular  to  it.  This  will  be  an  ellipse.  Then  the  two  wave- 
fronts  normal  to  OP  that  may  be  propagated  through  the  crystal  have  re¬ 
fractive  indices  equal  to  the  semi -axes,  OA  and  OB,  of  this  ellipse. 

The  displacement  vector  D  in  the  plane  polarized  wave  with  refractive 
index  equal  to  OA  vibrates  parallel  to  OA.  Similarly,  the  displacement 
vector  in  the  wave  with  refractive  index  equal  to  OB  vibrates  parallel 
to  OB.  From  this  it  follows  that  the  two  possible  waves  with  wave  normal 
Xj  have  refractive  indices  n^  and  n^  and  D  is  parallel  to  x^  and  x^ 
respectively.  Similarly,  the  wave  normal  x^  corresponds  to  two  waves 
of  refractive  indices  n^  and  n^,  with  D  parallel  to  x^  and  x^.  A  similar 
statement  applies  to  x3.-  For  this  reason,  n^,  n£,  and  n^  are  called  the 
principal  refractive  indices. 

Note  that  the  refractive  index  itself  is  not  a  tensor,  although  its 
variation  with  direction  is  determined  by  the  dielectric  constant  which 
is  a  tensor. 

Using  the  indicatrix,  crystals  can  be  classified  into  three  categories 
according  to  the  relationship  among  the  three  principal  indices  of  re¬ 
fraction. 

1 .  *  n^  *  n3 

The  ellipsoid  becomes  a  sphere.  For  any  wave-normal  direction, 
the  central  section  is  a  circle.  Thus,  there  is  no  natural 
birefringence  for  any  wave  normal  direction.  This  is  the  case 
in  cubic  crystals  and  isotropic  materials. 
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i  nj 

Here  only  one  section  Is  circular.  So,  for  only  wave-normal 
direction  Is  there  no  birefringence.  The  axis  Is  the  direction 
of  the  unique  index  of  refraction  is  called  the  optic  or  princi¬ 
pal  axis.  This  type  of  crystal  is  known  as  uniaxial.  Hexagonal, 
trigonal,  and  tetragonal  crystals  belong  to  this  group  as  do 
nematic  liquid  crystals. 

nl  *  n2  *  n3 

Now,  there  are  two  different  optic  axes.  The  crystal  is  called 
biaxial.  The  three  remaining  crystal  systems;  orthorhbmbic, 
monoclinic,  and  triclinic  belong  to  this  class. 
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APPENDIX  4C 


Nematic  Director  Orientation  To 
Match  Index  of  Refraction 
Variation 

Liquid  crystals  are  characterized  by  an  orientational  order 
of  their  constituent  rod-like  molecules  In  nematic 

liquid  cyrstals  this  orientational  order  has  uniaxial  symmetry, 
the  axis  of  uniaxial  symmetry  being  parallel  to  a  unit  vector, 
n,  called  the  director.  The  indicatrix  formulation  given  in 
Appendix  4B  applies  but  since  the  orientation  of  the  liquid 
crystal  molecules  can  be  altered  under  application  of  an 
external  field  the  orientation  of  the  indicatrix  can  be  changed. 
This  is  analogous  to  reorienting  a  solid  so  that  light  passing 
through  it  "sees"  a  different  index  of  refraction.  Figure  4.2-2 
illustrates  this  phenomenon.  In  this  example,  a  liquid  crystal 
possessing  positive  dielectric  anisotropy  is  shown  at  the  top 
of  the  figure  in  the  homogeneous  state.  Correspondingly,  the 
indicatrix  is  aligned  so  that  a  wave  polarized  in  the  x-direction 
is  influenced  by  an  index  of  refraction  ng.  At  the  bottom  of 
the  figure,  the  orientation  has  been  altered  to  the  homeotropic 
state.  In  which  case,  the  same  polarized  light  encounters  an 
index  of  refraction  nQ.  At  an  intermediate  point,  the  index 
of  refraction  is  dependent  on  9. 

4C.1  The  Value  of  the  Director 

In  Appendix  4A,  the  required  index  of  refraction  modulation 
to  match  an  ideal  thin  lens  was  derived.  To  create  that  modula¬ 
tion,  it  is  necessary  to  alter  the  orientation  of  the  liquid 
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crystal  molecules. 

Consider  the  slice  of  the  indicatrix  shown  in  Figure  AC-  1 

If  light  was  propagating  in  the  negative  yn  direction  and  polarized 

in  the  x  direction,  it  would  "see"  an  index  of  refraction  n  . 
n  e 

If  the  indicatrix  was  rotated  through  an  angle  Gn,  That  sane  ray 
would  "see"  an  index  of  refraction,  rn. 

The  slice  of  the  indicatrix  shown  can  be  expressed  as 


Equivalently, 


From  the  geometry  of  Figure  4C-J- 


(4C-1) 


(AC-2) 


(AC-3) 


(AC-A) 


(AC-5) 


(4C-6) 
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Or  finally, 


(4C-7) 


For  any  desired  index  of  refraction  between  nQ  and  ng,  Equation 

4C-7  gives  the  required  6  or  orientation  of  the  nematic  crystal 

director.  Note  that  indeed  when  r  *  n  ,  6  =  0°  and  when 

n  e 

r  =  n  ,  P  =  90°. 
n  o 

Now  the  results  of  Appendix  4A  can  be  factored  into  the  re¬ 
quired  orientation.  From  Appendix  4A, 

2 

r 

r  «=  n  - 

n  e  2Af  (4C-8) 

Where  terms  have  been  altered  to  match  their  meaning  here. 

The  required  form  for  the  index  of  refraction  can  now  be 
substituted  into  equation  (4C-7)  yielding  the  required  angle 
variation  for  the  liquid  crystal  director  as  a  function  of  the 
radial  distance  from  the  center  of  the  liquid  crystal  lens. 

As  an  example,  let, 
ne  -  1.634 

nQ  =  1.492 

A  *  14  x  lo'6  meters 
f  -  10  meters 

Figure  4C- 2,  shows  the  required  director  variation. 

Now  that  the  required  orientation  of  the  director  has  been  found, 
the  next  step  is  to  identify  what  field  is  required  to  produce  this 
variation. 


REQUIRED  DIRECTOR  ANGLE  DEGREES 
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Figure  AC-2  Radial  Distance  cm  x  10”^ 
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APPENDIX  4D 

Orientation  of  the  Liquid  Crystal 
Director  to  an  Applied 
Electric  Field 

For  most  display  applications,  liquid  crystals  are  used  as 
two  state  devices.  That  is,  the  crystal  is  either  in  its  equili¬ 
brium  state  or  is  driven  by  application  of  a  field  to  a  highly 
altered  state.  To  accomplish  this  result,  devices  are  normally 
operated  at  two  to  five  times  the  threshold  voltages.  Thus 
insuring  the  desired  effect. 

In  attempting  to  create  a  lens  with  a  liquid  crystal,  the 
transition  state  or  partial  reaction  of  the  crystal  to  an  applied 
field  is  critical.  This  appendix  described  the  orientation  of 
local  directors  to  an  applied  field. 

The  intial  formulation  follows  Sheng  (Chapter  8,  P.ef.11), 
Consider  a  very  small  spatial  region  inside  a  macroscopic  sample 
of  nematic  liquid  crystal  that  contains  a  sufficiently  large 
numbers  of  molecules  of  the  same  orientation  eo  that  a  local 
director  can  be  used  to  describe  this  orientation.  Imagine 
the  division  of  a  macroscopic  sample  into  such  small  spatial 
regions.  In  each  region,  define  an  orientational  director,  n, 
so  that  a  director  field  n(r)  results. 

Now  the  total  free-energy  density  of  the  system  under  external 
fields  can  be  written 

F  "  Fd  +  Fm  +  Fe  <*D-D 

Where 

Fj  -  Free-energy  density  associated  with  distortion  from 
equilibrium 
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Fm  -  Magnetic  field  contribution 
Fg  *  Electric  field  contribution 

The  total  free  energy  of  the  smaple  is  given  by 

F  -  /f  d3r  (4D-2) 

volume  of 
the  sample 

Sheng  expands  equations  (4D-1)  giving 

F  =  K1XC  v.n(r)]2  +  K22C|n(r)  .  VXn(r)  |  -  f  ]2 


+  K33C|n(r)  x  Vxn(r)  |  ]2  -  Ax  [8.n(r)  ]2  -  ^-AeCS  .n  (r)  ]2} 


(AD-3) 


Where 

Kn  «=  Spray  elastic  constant 
K22  "  Twist  elastic  constant 

K33  *  Bend  elastic  constant 

4*  =  Xu  *  Xj_ 

A  ^ 

Xii  ■  Susceptibility  per  unit  volume  parallel  to  n(r) 

X  ■  Susceptibility  per  unit  volume  perpendicular  to  n(r) 

Ae  -[en  -eje0 

e,,  *  Dielectric  constant  parallel  to  n(r) 

A 

Ej^  ■  Dielectric  constant  perpendicular  to  n(r) 
eQ  -  Permittivity  of  free  space 

H  *  Externally  applied  magnetic  field 

E  -  Externally  applied  electric  field 
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Figure  4D-1  gives  an  illustration  to  the  meaning  of  the  first  three 
terms  on  the  right  hand  side  of  Equations  (4D-3) . 

By  knowing  that  the  liquid  crystal  molecules  will  orient  to  mini¬ 
mize  the  free  energy  (Equation  4D-2),it  is  now  possible  to  use 
Equation  (4D-3)  and  find  the  orientation  of  the  liquid  crystal 
director  as  a  function  of  applied  field  strength  and  position 
within  the  sample. 

Significant  mathematical  simplification  results  if  the  three 
constants,  K^,  K22'  an<*  K33'  co^^ective^y  called  the  Frank  con¬ 
stants  are  set  equal  to  K.  Similarly  if  the  -it / x q  term  is  omitted 
which  does  not  add  to  the  minimization  problem,  Equation  (4D-3) 
can  be  simplified  to 

F  -  |  { [ V .n(r)  3 2  +  C|v  x  n(?)|l2  -  ^tS.J(r)]2  -  ^j-^CE.nC?)  3  > 

(4D-4) 

Consider  the  geometry  in  Figure  (4D-2)  .  A  positive  dielectric 

liquid  crystal  fills  a  space  of  width  d  centered  at  x  =  0.  Molecules 

are  constrained  at  the  surface  to  be  parallel  to  the  enclosure 

walls.  With  the  application  of  an  external  electric  field  in  the 
x  direction,  molecules  away  from  the  walls  re-orient  and  attempt 

to  line  up  with  the  field.  From  Figure  (4D-?)  the  components 
of  the  local  director,  can  be  written  as, 

nx  -  cosG(x) 

n y  ■  0 


sinO(x) 


(4DO) 
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Figure  4D-1  Three  types  of  distortion  in  a  director  field: 

(a)  gives  V.ft(r)^  0,  (b)  gives  ft(r).Vx£(r)  ^  0, 
and  (c)  gives  n(r)  x  Vxn(r)  f  0.  Each  director 
is  shown  with  double  arrows  in  order  to  indicate 
that  the  "head"  and  the  "tall"  directions  represent 
exactly  the  sane  physical  state  in  a  nematic 
sample. 
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Figure  4D-2  Molecule  Orientation  With  Applied  Field 
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From  which 


A  O 

CV.n(r) j 


sin2© 


(4D-6) 


[  |Vxn(r)  |  ]2  =  cos2© 


(4D-7) 


And 


And 


2^  AeCE.n(r) ]2 


"Ay  CH.n(r)]2  -  0 


-E  ie  2~ 
=2~~  cos  0 


(4D-8) 


(4D-9) 


Substitution  of  these  into  equation  (4D-4)  gives  the  free 
energy  density  of  the  system  as 


v  „  K,rd©  2  E2Ae  „  2a, 

F  "  7{CH£3  '  OT-  cos  0} 


Correspondingly,  the  total  free  energy  is  given  by 
v  m  t  K,rd©  2  *2h£  ooc 2a 

r  J  -  /. cos  0  ;dr 


(4D-10) 


4ttK 


(4D-11) 


It  is  this  function  that  the  molecules  move  to  minimize.  To 
minimize  an  integral  of  the  form 
b 


/  G(y(x),  **£>■,»)  dx 


(4D-12) 


By  varying  y(x),  the  optimal  function  y(x)  must  satisfy  the 
equation 


dG  d  r  3G 


•>  0 


(4D-13) 


V . 


k&&'EBSi3S>i 


Which  is  called  the  Euler-Lagrange  Equation.  Comparing  Equations 
(4D-11)  and  (4D-12)  and  noting  that  the  only  variations  in  the 
volume  is  in  the  x-direction  gives 


G  “Cal:i  "  hvr  cos2g 

And 

y  -  6 

Therefore,  the  Euler-Lagrange  Equation  becomes 

Ifir1  cos0sine  -  12  ai ]  =  0 

Equivalently, 

p.2  j2~ 

cosQsine  -  -  0 


For  simplicity,  let 


So 


(4D-14) 


(4D-15) 


(4D-16) 


(4D-17) 


cosAsin© 


(4D-18) 
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Rewriting 


2 

sinOcosO  d©  _  d  9 •  d©  _  1  d  rd©\ 
5"7  3x  dx  dx  -  7  dx  ',dx' 

E 


(4D-19) 


sin©  dsin©  1  d  /C19n2 

JT~  ST  =  7 

^  E 

Or 

1  1  dsin2©  1  d  ,d©N2 

7  f7  “dx  .  =  7  dx  W 

S  E 

Integration  Yields 


(4D-20) 


(4D-21) 


<d©) 

die 


=  (sin©)  + 

%E 


(4D-22) 


The  constant  of  integration  can  be  evaluated  by  noting  that  be¬ 
cause  of  symmetry  4^1  =0.  Calling  the  corresponding  0(x=O)=0 

ax  x=0  ’ 


the  minimum  value  of  0 


C  =  - 


,  sin  Qm 

(  - - 

b  E 


0 


(4D-23) 


So  equation  (4D-22)  becomes 

(dg  )  2  .(  sln£  ,2  _(3inem)2  (4D-24) 

-5- 

Equivalently, 

f  -  t  -  sln\ 


(4D-25) 
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Where  the  negative  sign  applies  for  -d/24x4  0  and  the  positive 
sign  applies  for  0  <  x  4  d/2.  Because  of  symmetry  only  one 
section  needs  to  be  considered.  So,  for  -d/24  x4  0, 

e 

-  sin2em  (4D-26) 


de 

Sinze'  -  sin^e^  (4D-27) 


Equivalently, 


0 


/2 


Unfortunately,  there  is  no  closed  form  solution  to  the  right  side. 
Insight  can  be  gained,  however,  by  assuming  that  sufficient  field 
strength  has  been  applied  to  drive  0m  to  zero.  In  that  case. 
Equation (4D-27)  becomes 


fE  lK  +  f]  ’  •  / 
»/2 


de  , 

sin^G' 


5 

v/2 


dQ 

s  inG 


Or 


—  jx  +  jJ  m  log  tan  j 


fi 

Finally 

0-2  Tan 


1  expK 


(x  +  |) 


■d/2  <  x  <  0 


(4D-28) 


(4D-29) 


(4D-30) 
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As  an  example  let 


d 

K 

A  E 


14  x  lo"  meters 

10“6  dynes  =  10”^  newtons 

.7  x  .  -  6.195  x  10-12 


Figure4D-*3  shows  the  results  for  two  different  voltages  across  the 
cell.  One  very  interesting  fact  comes  out  of  this  example.  The 
average  orientation  of  the  directors  are 

ff50  -  22.5° 

And 

5100  "  14*9° 

If  n_  =  1.64  and  n  -  1.492,  then  the  average  indices  of  refraction 

c  O 

are 


This  is  a  small  change  in  index  of  refraction  for  a  significant 
change  in  orientation  and  a  very  significant  change  in  applied  voltage. 


DIRECTOR  ORIENTATION  0E6REES 
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Figure  4D-3  Director  Orientation  Across  Liquid  Crystal 
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